Effective field theory for models defined over small-world networks. 
First and second order phase transitions. 
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We present an effective field theory method to analyze, in a very general way, models defined over 
small-world networks. Even if the exactness of the method is limited to the paramagnetic regions 
and to some special limits, it provides, yielding a clear and immediate (also in terms of calculation) 
physical insight, the exact critical behavior and the exact critical surfaces and percolation thresholds. 
The underlying structure of the non random part of the model, i.e., the set of spins filling up a given 
lattice Co of dimension do and interacting through a fixed coupling Jq, is exactly taken into account. 
When Jo > 0, the small-world effect gives rise, as is known, to a second-order phase transition that 
takes place independently of the dimension do and of the added random connectivity c. When 
Jo < 0, a different and novel scenario emerges in which, besides a spin glass transition, multiple 
first- and second-order phase transitions may taJte place. As immediate analytical applications we 
analyze the Viana-Bray model (do = 0), the one dimensional chain (do = 1), and the spherical 
model for arbitrary do. 

PACS numbers; 05.50.-|-q, 64.60.aq, 64.70.-p, 64.70.P- 



I. INTRODUCTION 

Since the very beginning of the pioneer work by Watts 
and Strogatz the interest toward small- world net- 
works - an interplay between random and regular net- 
works - has been growing "exponentially" . Mainly, there 
are two reasons that have caused such a "diffusion" . 

The first reason is due to the topological properties of 
the small- world network. In synthesis, if N is the size 
of the system, for any finite probability p of rewiring, or 
for any finite added random connectivity c (the two cases 
correspond to slightly different procedures for building a 
small-world network) one has: a "short-distance behav- 
ior" , implying that the average shortest distance between 
two arbitrarily chosen sites grows as 1{N) ~ log(A^), as 
in random networks, and a large clustering coefficient, 
C{N) ~ 0(1), as in regular lattices. The interplay be- 
tween these two features makes small- world networks rep- 
resentative of many realistic situations including social 
networks, communications networks, chemical reactions 
networks, protein networks, neuronal networks, etc. 

The second reason is due the fact that, in models de- 
fined over small-world networks, despite the presence of 
an underlying finite dimensional structure - a lattice Cq 
of dimension do < oo - the existence of short-cut bonds 
makes such models mean-field like and - hopefully - ex- 
actly solvable. However, even if such a claim sounds in- 
tuitively correct, the complexity of these models turns 
out to be in general quite high and, compared to numeri- 
cal works, there are still few exact results for small- world 
networks 0-H (for t he p ercolation and synchronization 
problem see [Ifll and pl|). 

In particular, for do > 1 a mean-field critical behav- 
ior is expected and has been also supported by Monte 
Carlo (MC) simulations Some natural questions 

then arise. Are we able to prove analytically such a be- 



havior? If for example do — 2, does the mean- field critical 
behavior hold for any situation? Yet, does the correlation 
length diverge at the critical temperature? 

Furthermore, even if for dg = 1 an exact analytical 
treatment has been developed at the level of replica sym- 
metry breaking (RSB) Q and one step replica symmetry 
breaking (IRSB) the calculations are quite involved 
and the solutions of the coupled equations to evaluate 
the order parameters require a certain numerical work, 
which becomes rapidly hard in the IRSB case. In any 
case, even if these methods are able to give in principle 
exact results at any temperature, they are not in general 
suitable to provide a clear simple and immediate physical 
picture of the model, even within some approximations. 
The main problem in fact resides in the presence of short 
loops: as soon as do > 1 these loops cannot be neglected 
and the "traditional" cavity and replica methods seem 
hardly applicable. In particular, we are not able to pre- 
dict what happens, for example, if we set Jq negative. 
Should we still expect a second-order phase transition? 
And what about the phase diagram? 

In this paper, we present a general method to study 
random Ising models defined on small-world graphs built 
up by adding a random connectivity c over an underlying 
arbitrary lattice Co having dimension do. We will show 
that this method - in a very simple and physically sound 
way - provides an answer to the above questions as well 
as to many others. 

Roughly speaking, as an effective field theory the 
method generalizes the Curie- Weiss mean-field equation 
m — tanh(/? Jto) to take into account the presence of the 
short-range couplings Jo besides the long-range ones J. 
As we will show, the magnetization m of the model de- 
fined over the small-world network, shortly the random 
model, behaves as the magnetization mg of the model 
defined over Co, shortly the unperturbed model, but im- 
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mersed in an effective external field to be determined 
self-consistently. Even if the exactness of this method is 
limited to the paramagnetic regions (P), it provides the 
exact critical behavior and the exact critical surfaces, as 
well as simple qualitative good estimates of the corre- 
lation functions in the ferromagnetic (F) and spin glass 
regions (SG). Furthermore, in unfrustrated systems, the 
method becomes exact at any temperature in the two 
limits c ^ 0+ and c — > oo. 

The consequences of such a general result are remark- 
able from both the theoretical and the practical point of 
view. Once the explicit form of the magnetization of the 
unperturbed model, uiq — mo{pjQ, /3h), is known, an- 
alytically or numerically, as a function of the couplings 
Jo and of the external field h, we get an approximation 
to the full solution of the random model, which is ana- 
lytical or numerical, respectively, and becomes exact in 
the P region. If we do not have mo = mo(/3Jo, but 
we know at least some of its properties, we can still use 
these properties to derive certain exact relations and the 
critical behavior of the random model. 

In the first part of the paper, after presenting the self- 
consistent equations, we focus on their application for 
a general study of the critical surfaces and of the criti- 
cal behavior. In the second part, we apply the method 
to study models of interest which can be solved analyti- 
cally (and very easily) as for them we know mo{/3JQ, ph): 
the Viana-Bray model (which can be seen as a do = di- 
mensional small- world model) , the one-dimensional chain 
small-world model, and the spherical small-world model 
in arbitrary do dimension. 

We stress that the critical surfaces as well as the corre- 
lation functions in the P region provided by the present 
method are exact and not based on any special ansatz as 
the replica-symmetry and the tree-like ansatz. We prove 
in particular that: independently of the added random 
connectivity c, of the underlying dimension da, of the 
structure of the underlying lattice Cq, and of the na- 
ture of the phase transition present in the unperturbed 
model (if any) , for Jq > we always have a second-order 
phase transition with the classical mean-field critical in- 
dices but with a finite correlation length if calculated 
along the "Euclidean distance" defined in Cq; on the 
other hand, for Jq < we show that, as soon as c is 
sufficiently large, there exist at least two critical temper- 
atures which, depending on the behavior of xoiPJoj Ph) 
- the susceptibility of the unperturbed system - corre- 
spond to first- or second-order phase transitions. This 
phenomenon will be explicitly shown in the example of 
the one-dimensional small-world model. Note that, as it 
will result from the detailed analysis of the self-consistent 
equation (Sec. IIIB), in any case, the critical behavior of 
the unperturbed model - if any - can never influence the 
behavior of the random model. 

The paper is organized as follows. In Sec. II we in- 
troduce the class of small- world networks over which we 
define the random Ising models, stressing some impor- 
tant differences concerning the definition of the correla- 



tion functions with respect to those usually considered 
in "ordinary" random models. In Sec. Ill we present 
our method: in Sec. IIIA we provide the self-consistent 
equations and their relations with physical correlation 
functions, in Sec. IIIB we analyze the stability of the 
solutions of the self consistent equations and the critical 
surface and behavior of the system. We separate the Sec. 
IIIB in the sub-cases Jo > and Jo < 0. In Sec. IIIC we 
discuss the limits of the method. In Sec. HID we study 
the stability between the F and the SG phases and the 
phase diagram. Finally, in Sec. HIE we mention how to 
generalize the method to cases with more different short- 
range couplings Jo, and to analyze possible disordered 
antiferromagnetic systems. In Sees. IV, V and VI the 
theory is applied to the three above mentioned example 
cases. The successive Sees. VII, VIII and IX are devoted 
to the derivation of the method. The starting point of 
the proof is given in Sec. VII and is based on a general 
mapp ing between a random model and a non random one 
suitably adapted to the present case. The self- 
consistent equations are then easily derived in Sec. VIII. 
Note that, apart from the equations concerning the sta- 
bility between the P-F and the P-SG transitions, which 
are derived in Sec. IX, the derivations of the equations 
presented in Sec. IIIB are mostly left to the reader, since 
they can be easily obtained by standard arguments of 
statistical mechanics using the Landau free energy ^{m) 
that we provide and that is derived in Sec. VIII too. Fi- 
nally, in Sec. X we draw some conclusions. In Appendix 
A we generalize the method to inhomogeneous external 
fields to make clear the subtle behavior of the correlation 
functions in small-world models. 



II. RANDOM ISING MODELS ON 
SMALL- WORLD NETWORKS 

The family of models we shall consider are random 
Ising models constructed by super-imposing random 
graphs with finite average connectivity c onto some given 
lattice Co whose set of bonds (i,j) and dimension will 
be indicated by Fo and do, respectively. Given an Ising 
model {the unperturbed model) of N spins coupled over 
Co through a coupling Jo with Hamiltonian 

Ho - Jo ^ (JiCTj -h'^cTi, (1) 

(»j)ero * 

and given an ensemble C of unconstrained random graphs 
c, c G C, whose bonds are determined by the adjacency 
matrix elements Ci_j = 0, 1, we define the corresponding 
small-world model as described by the following Hamil- 
tonian 

Hc,j'^=^ Ho -'^c.tjJij(7iaj, (2) 
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the free energy F and the averages (O)', with Z = 1,2, 
being defined in the usual (quenched) way as (/3 = 



and 



i^Oy'ilY.P^c) I dV{J){OY,^j, 1 = 1,2 (4) 

where Z^.j is the partition function of the quenched sys- 
tem 



Zo.j - J2 



-/3ffo,j({<T.}) 



(5) 



{0)c,j the Boltzmann-average of the quenched system 
((O) depends on the given reahzation of the J's and of 
c: (C) = (0)c;j; for shortness we later will omit to write 
these dependencies) 



(6) 



and dV (J) and P(c) are two product measures given 
in terms of two normalized measures dfi{Jij) > and 
p{c-i,j) ^ Oi respectively: 

dvij)"^ n y'dA'(j.,,)=i, (7) 



p{c)'= n ^'(^^^■), E Mc.,)=i. (8) 



The variables Cij € {0, 1} specify whether a "long- 
range" bond between the sites i and j is present (q j = 1) 
or absent (cij = 0), whereas the Jij 's are the random 
variables of the given bond (« , j) . For the Jij 's we will not 
assume any particular distribution, while, to be specific, 
for the Cij 's we shall consider the following distribution 



N 



(9) 



This choice leads in the thermodynamic limit — > cxd to 
a number of long range connections per site distributed 
according to a Poisson law with mean c > (so that in 
average there are in total cN/2 bonds). Note however 
that the main results we report in the next section are 
easily generalizable to any case in which Eq. ^ holds, or 
holds only in the thermodynamic limit due a sufhcicntly 
small number of constrains among the matrix elements 

CiJ. 

When we will need to be specific, for the Jij 's we will 
assume either the distribution 



dJ, 



5 {Ji^j - J) , 



(10) 



dfj.{J^j) 



HF =^ J2 ^(^) / (^-■■^) (3) dJ.^,J 

cec 



p6 {.h,j - J) dJ,,, + (1 - p)S ( + J) , (11) 



to consider ferromagnetism or glassy phases, respectively. 
InEq. [0,1]. 

The quantities of major interest are the averages, and 
the quadratic averages, of the correlation functions which 
for shortness will be indicated by C^^^ and C^^-' . For ex- 
ample, the following are non connected correlation func- 
tions of order k: 



(12) 
(13) 



where fc > 1 and the indices ii, . . . ,ik are supposed all dif- 
ferent. For shortness we will keep on to use the symbols 
C*-^-* and C*-^^ also for the connected correlation function 
since, as we shall see in the next section, they obey to 
the same rules of transformations. We point out that 
the set of indices ii, . . . ,ik is fixed along the process of 
the two averages. This implies in particular that, if we 
consider the spin with index i and the spin with index j, 
their distance remains undefined, or more precisely, the 
only meaningful distance between i and j, is the distance 
defined over Cq, i.e., the Euclidean distance between i 
and j, which we will indicate as ||i — j Hq. 

Therefore, throughout this paper, it must be kept in 
mind that, for example, C*^"'^^(||i — j|||^) = {cicrj) is very 
different from the correlation function G^^-'(Z) of two 
points at a fixed distance I, I being here the distance 
defined over both and the random graph c, i.e., the 
minimum number of bonds to join two points among both 
the bonds of Fq and the bonds of the random graph c. 
In fact, if, for Jq = 0, one considers all the possible re- 
alizations of the Poisson graph, and then all the possible 
distances I between two given points i and j, one has 



N 

= J2Pn{1)G^'\1) (14) 



1=1 



where here Pn{1) is the probability that, in the system 
with N spins, the shortest path between the vertices i 
and j has length I. If we now use G''^'(Z) ~ (tanh(/3J))' 
[l6l | (in the P region holds the equality) and the fact that 
the average of I with respect to Pn{1) is of the order 
log(iV), we see that the two point connected correlation 
function (|14p goes to in the thermodynamic limit. Sim- 
ilarly, all the connected correlation functions defined in 
this way are zero in this limit. Note however, that this 
independence of the variables holds only if Jq = 0. This 
discussion will be more deeply analyzed along the proof 
by using another point of view, based on mapping the 
random model to a suitable fully connected model. 
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III. AN EFFECTIVE FIELD THEORY 
A. The self-consistent equations 

Depending on the temperature T, and on the param- 
eters of the probabiUty distributions, dfi{-) and p{-), the 
random model may stably stay either in the paramag- 
netic (P), in the ferromagnetic (F), or in the spin-glass 
(SG) phase. In our approach for the F and SG phases 
there are two natural order parameters that will be indi- 
cated by m'-^-' and m^^'^\ Similarly, for any correlation 
function, quadratic or not, there are two natural quanti- 
ties indicated by C^^' and C^^*^', and that in turn will 
be calculated in terms of m^^^ and m'^^'^\ respectively. 
To avoid confusion, it should be kept in mind that in our 
approach, for any observable O, there are - in principle 
- always two solutions that we label as F and SG, but, 
as we shall discuss in Sec. HID, for any temperature, 
only one of the two solutions is stable and useful in the 
thermodynamic limit. 

In the following, we will use the label o to specify that 
we are referring to the unperturbed model with Hamil- 
tonian ([T]). Note that all the equations presented in this 
paper have meaning and usefulness also for sufficiently 
large but finite size N. For shortness we shall omit to 
write the dependence on N. 

Let mQ[(3JQ, ph) be the stable magnetization of the 
unperturbed model with coupling Jq and in the presence 
of a uniform external field h at inverse temperature (3. 
Then, the order parameters m'^^\ S=F,SG, satisfy the 
following self-consistent decoupled equations 



m 



) =TOo(/3J^^\/3j(5^)m(^) (15) 



where the effective couplings J^^\ j'^*^', and Jq^'^'' 
are given by 



/3j(F) / d^(J,j)tanh(/3J,j), 



/3J(SG) / d^(j;^j.)tanh^(/3J,j), 



and 



t(F) _ J 



/3jf = tanh"^(tanh^(/3Jo)). 



(16) 



(17) 



(18) 



(19) 



Note that > 4^'^\ 



For the correlation functions C'^^\ E=F,SG, for suffi- 
ciently large N we have 

C(^) = Co(/3jf \/3J(^)m(^) -^/g/i) -f O (^1) , (20) 

where Cq (/? Jq i Z^^) is the correlation function of the un- 
perturbed (non random) model. 



Concerning the free energy density / we have 
J d/i(J,,,)log[cosh(/3J,,,)] 
^Ihn^^ 5Z log[cosh(/3Jo)]-log[2cosh(/3/i)] 



N 
1 



+{J™ M E cosh(/3jf)) 
('■j)ero 

+ log[2cosh(/3/i)]} X i + iL(s)(m(^)), 



(21) 



where / = 1, 2 for I]=F,SG, respectively, and the non 
trivial free energy term is given by 

foi/3Jo, (3h) being the free energy density in the thermo- 
dynamic limit of the unperturbed model with coupling 
Jo and in the presence of an external field h, at inverse 
temperature p. 

For given f3, among all the possible solutions of Eqs. 
(|15p . in the thermodynamic limit, for both I]=F and 
S=SG, the true solution m^^\ or leading solution, is the 
one that minimizes L^^^: 

L^^'fm^^)) = min L^^) (m) . (23) 

V / me[-l,l] 

Finally, let k be the order of a given correlation func- 
tion C*^^-* or C*-^^. The averages and the quadratic aver- 
ages over the disorder, C^^^ and C^^^, are related to C^^^ 
and C*^^*^-* , as follows 

Cd) ^ C7(F), inF, (24) 
C^i) = 0, /fcodd, inSG, (25) 
^(1) ^ (j{SG)^ keven, in SG, (26) 



and 



C(2) 



, inF, 
C^^^'V, inSG. 



(27) 
(28) 



From Eqs. (H?]) and (gSl) for fc = 1, we note that the 
Edward- Anderson order parameter [l3| C*^^-* — (c)^ — 
QEA is equal to {C'^^^'')'^ = {m^^^^ f only in the SG 
phase, whereas in the F phase we have qea = (^tt-*-^-*)^. 
Therefore, since m^^^) ^ m^^\ m^^*^^ is not equal to 
y/QEA', in our approach m*^^'^-' represents a sort of a spin 
glass order parameter [l^ . 

The localization and the reciprocal stability between 
the F and SG phases will be discussed in Sec. HID. 
Note however that, at least for lattices Cq having only 
loops of even length, the stable P region is always that 
corresponding to a P-F phase diagram, so that in the P 
region the correlation functions must be calculated only 
through Eqs. ^ and ((^ . 
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As an immediate consequence of Eq. dTSl) we get the 
susceptibility x^^^ of the random modeh 



Xo 



(3h 



PJ^^^Xo (^l3jP,l3J(^)mi^) +I3h 



1 



, (29) 



where xo stands for the susceptibihty xo of the unper- 
turbed model divided by /3 (we will adopt throughout 
this dimensionless definition of the susceptibility) : 

~ ^flT oi^s def dmo{pJo,f3h) 1 dniQ {pjo, f3h) 

^ dm ^ -p — dh — 

and similarly for the random model. 

For the case S =:F without disorder {dfi(J') — d{J' — 
J)dJ'), Eq. ([29]) was already derived in Q by series 
expansion techniques at zero field {h = 0) in the P region 
(where to = 0). 

Another remarkable consequence of our theory comes 
from Eq. I|20p. We see in fact that in the thermodynamic 
limit any correlation function of the random model fits 
with the correlation function of the unperturbed model 
but immersed in an effective field that is exactly zero in 
the P region and zero external field {h = 0). In other 
words, in terms of correlation functions, in the P region, 
the random model and the unperturbed model are in- 
distinguishable (modulo the transformation Jq Jq 
for E =SG). Note however that this assertion holds only 
for a given correlation function calculated in the ther- 
modynamic limit. In fact, the corrective 0(l/iV) term 
appearing in the rhs of Eq. ([^0]) cannot be neglected 
when we sum the correlation functions over all the sites 
I € £o, as to calculate the susceptibility; yet it is just this 
corrective 0(1/A^) term that gives rise to the singularities 
in the random model. 

More precisely, for the two point connected correlation 
function 



where Z = 1, 2 for E = F, SG, respectively, if 



(31) 



def I , 



(32) 



we have 



x!? = Xo..(/3jf\/3/^)™(^^+/3M + ^ 



(33) 



where the dependence on TV in xl j and xo are under- 
stood. Eq. pS]) clarifies the structure of the correlation 
functions in small-world models. In the rhs we have two 
terms: the former is a distance-dependent short-range 
term whose finite correlation length, for T ^ T^q (T^q 



being the critical temperature of the unperturbed model 
with coupling Jg ), makes it normalizable, the latter 
is instead a distance-independent long-range term which 
turns out to be normalizable thanks to the 1/A^ factor. 
Once summed, both the terms give a finite contribution 
to the susceptibility. It is immediate to verify that by 

summing xl / over all the indices i,j G Co and dividing 
by TV we get back - as it must be - Eq. Eq. ((55)) 

will be derived in Appendix A where we generalize the 
theory to a non homogeneous external field. 



B. Stability: critical surfaces and critical behavior 



Note that, for /? sufficiently small (see later), Eq. pS]) 
has always the solution to*^^-' = 0, and furthermore, if 
TO^^^ is a solution, — to^^^ is a solution as well. From 
now on, if not explicitly said, we will refer only to the 
positive (possibly zero) solution, the negative one being 
understood. A solution to*^^^ of Eq. ([T3]) is stable (but 
in general not unique) if 



1 - 



/SJ^^^Xo (/3J(i^\/3j(^^TO(^) +/3/i) > 0. (34) 



For what follows, we need to rewrite the non trivial 
part of the free energy density L^^^ (m) as 



(m) = f3fo (/? Jf ^ , O) - TOO (/? Jf ^ , O) Ph 

+iP^^'^ (to) , (35) 



where the introduced term -i/j^^^ plays the role of a Lan- 
dau free energy density and is responsible for the critical 
behavior of the system. Around to = 0, up to terms 
0{h^) and Oim^h), ip'''^\m) can be expanded as follows 



V'(^) (to) 



2 4 6 



where 



1 - 



+A(/3/o)(/3jf' ,/3J(^)to), 
/3J(^'X0 (/3jf\0)]/3j(^), 



= 



d{ph)- 



diphy 



■Xo 



■Xo 



3! 



5! 



(36) 



(37) 



(38) 



(39) 



/i(^) = mo \o) J(^) +X0 \o) /3/^)/3/i(40) 
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finally, the last term A (/3/o) (/?Jo, /3j'^^m) is defined 
implicitly to render Eqs. ([55]) and ([55]) exact, but terms 
Oih?) and 0{m^h), explicitly: 



A 



(/5/o) (/3jf\/5j(^'r 



E 



32fc-2 



2fe-2 



2fe 



(2fc)! 



-{.41) 



We recall that the fc — 2-th derivative of xo {pj^^^ ■, 
with respect to the second argument, calculated at 
h = 0, gives the total sum of all the A:-th cumu- 

lants normalized to N: d^^"^ xq {j3 (3hj \h=o = 

T.^^,...,^S^il■■■(^ik)Q^ IN, where (cr^i • • • cr^Jo^^ stands 
for the cumulant, or connected correlation function, of or- 
der k of the unperturbed model, {cfiiCfi.^)^^^ = {oi^ai^)^) — 
(o'ii)o(o'i2)o, etc.. Note that, apart from the sign, these 
terms are proportional to the Binder cumulants [Toj 
(which are all zero above Tco for fc > 2) only for N finite. 
In the thermodynamic limit the terms b'^^\ . . ., in 

general are non zero and take into account the large de- 
viations of the block-spin distribution functions from the 
Gaussian distribution. 

Let Tc = l//3c be the critical temperatures, if any, 
of the random model and let t^^^^ be the corresponding 
reduced temperatures: 



(s) T — Tc 



/3f ) - /? 



T, 



13, 



o{t 



(S)^,2 



(42) 



Here, the term "critical temperature" , stands for any 
temperature where some singularity shows up. However, 
if we limit ourselves to consider only the critical temper- 
atures crossing which the system passes from a P region 
to a non P region, from Eq. (l34l) it is easy to see that, 
independently on the sign of Jq and on the nature of the 
phase transition, we have the important inequalities 



/3F 



^ l-'cO ' 



(43) 



where we have introduced (3^q , the inverse critical tem- 

perature of the unperturbed model with coupling Jq 
and zero external field. If more than one critical tem- 
perature is present in the unperturbed model, /3^q is the 
value corresponding to the smallest value of these crit- 
ical temperatures (highest in terms of /?). Formally we 
set /3^Q = cx) if no phase transition is present in the un- 
perturbed model. A consequence of Eq. is that, in 
studying the critical behavior of the system for h = 0, 
we can put /i'^^ — 0. Throughout this paper, we shall 
reserve the name critical temperature of the unperturbed 
model as a P-F critical temperature through which the 
magnetization toq {f3Jo, 0) passes from a zero to a non 
zero value, continuously or not. This implies, in particu- 
lar, that for Jo < we have - formally - (3co — oo. 

In this paper we shall study only the order parameters 
m^^^ and m^^*^-* , whereas we will give only few remarks on 



how to generalize the method for possible antiferromag- 
netic order parameters. We point out however that the 
existence of possible antiferromagnetic transitions of the 
unperturbed model does not affect the results we present 
in this paper. 

It is convenient to distinguish the cases Jo > and 
Jq < 0, since they give rise to two strictly different sce- 
narios. 



The case Jq > 



In this case PJ^^^Xa {f3Jo,PJ'-^^m'^^^+l3h) 



IS an m- 

creasing function of /3 for /3 < /3^q (and for h>0). As a 
consequence, we have that for sufficiently low tempera- 
tures, the solution m*^^^ = of Eq. becomes unsta- 
ble and two - and only two - non zero solutions ±m*-^^ are 
instead favored. The inverse critical temperatures (3, 



(F) 



(SG) 

and f3c can be determined by developing - ioi h = - 
Eqs. (fT5|l for small m'^^ and m^^'^\ respectively, which, 
in terms of xo gives the following exact equation 



Xo(/3f)jf\0)/3f)j(^) = l, /3F 



(44) 



where the constrain Pc < Pco excludes other possible 
spurious solutions that may appear when do > 2 (since 

in this case /3^q^ may be finite). 

The critical behavior of the system can be derived by 
developing Eqs. for small fields. Alternatively, one 
can study the critical behavior by analyzing the Landau 
free energy density tp^^^m^^^) given by Eq. ([36|) . 

In the following we will suppose that for Jo > 0, b^^' 
be positive. We have checked this hypothesis in all the 
models we have until now considered and that will be 
analyzed in Sees. IV, V and VI. Furthermore, even if the 
sign of c*-^-* cannot be in general a priori established, for 
the convexity of the function /o with respect to (3h, the 

sum of the six-th term with A (/3/o) (/? J($^^ , f3J^^^ m(^)) , 

in Eq. (|36p must go necessarily to +oo for m*-^) — > cx). In 
conclusion, when Jo > 0, for the critical behavior of the 
system, the only relevant parameters of V'*'^'' are a'-^-', 
6(^) and /i(^) = x{f3Jo^\0)J^^^h, so that the critical 
behavior can be immediately derived as in the Landau 
theory for the so called model [lO]. On noting that 



a(^) > 0, for i(^) > 0, 
a(^' < 0, for < 0, 



it is convenient to define 



op 



(45) 



(46) 



so that we have 



a(^)=A(^)| ,,j(^)-fO(t(^))^ (47) 

13= 
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Note that, due to the fact that Jo > 0, A^^) > 0, 
and, as already mentioned, fe^^-* > as well. By using 
Eq. (|47| for (3 < Pco and near (3c , we see that the 
minimum ?fi^^^ of V'^^^ the solution of Eq. (jl5p 
near the critical temperature, is given by 

[0, t(^) >0, 

Similarly, we can write general formulas for the sus- 
ceptibility and the equation of state. We have 



.(.,7IW+0(1), t(^) <0, 



(49) 



/l3 +0 



(/.i)(50) 



Finally, on using Eqs. and (|48p we get that the 

specific heat C^^-* has the following finite jump disconti- 
nuity at /3, 



where Cc is the continuous part of the specific heat cor- 
responding to the part of the free energy density without 

Hence, as a very general result, independently of the 
structure of the underlying graph and its dimension 
c?o, independently of the nature of the phase transition 
present in unperturbed model (if any), and independently 
of the added random connectivity c, provided positive, 
we recover that the random model has always a mean- 
field critical behavior with a second-order phase transi- 
tion with the classical exponents /3 — ^/l, 7 = 7' = 1, 
5 — Z and a = a' = 0, and certain constant coefficients 
depending on the susceptibility xo and its derivatives cal- 

culated at /3 = /3c and external field h — Q. Note how- 
ever, that the correlation length of the system calculated 
along the distance of Cq, \\ ■ ||o, remains finite also at 
fic ■ In fact, from Eq. (PO)) . for the two point correla- 

def 

tion function at distance r — \\i — in Co we have 

C(^)(r) =Co(/3jf\/3/^)m(^) -f/3/i;r). (52) 

If we now assume for Co(/3Jo, 0; r) the following general 
Ornstein-Zernike form 



t(^) >0, 
t(^) < 



(51) 



Co(/3Jo,0;r) 



(53) 



fo ('') = /o (PJo ; J') being a smooth function of r (which 
has not to be confused with the free energy density), and 



^0 = Co(/3Jo) the correlation length, which is supposed 
to diverge only at /3co (if any), on comparing Eqs. ([5^ 
and ([55]) for [3 > (ic we have (notice that, as explained 
in Sec. IIIA, at least for lattices Cq having only loops 
of even length, the physical correlation function is only 
that corresponding to E =F, i.e., C^^^ = C^^^) 



C(^)(r) = 



where 



and 



/^'^^(0-/o(rf^r-), 



Therefore, due to the inequalities p3|) . we see that 



(S)| 



l/3=/3i 



= ^o(/3f^jf^)<oo. 



(54) 



(55) 



(56) 



(57) 



The knowledge of Co(/3Jo, Ph; r) also for h ^ would al- 
lows us to find the general expression for C^^-* (r) through 

Eq. §^ also for P > (31'. However, since Co(/3Jo, I3h] r) 
has no critical behavior for /i ^ 0, it follows that C*^^^ (r) 
cannot have a critical behavior for (3 > (3c either (and 
then also for (3 ^ (3c from the right). This result is 
consistent with |2|]. 



In this case J', 



(F) 



The case Jo < 



< 



0, so that - in general - 
longer a monotonic func- 
tion of (3. However, it is easy to see that that for (3 = 
and (3 00, this function goes to 0. Therefore, for a 
sufficiently large connectivity c, from Eq. we see 

that there may appear at least two regions where the 
paramagnetic solution m^^^ = is stable, separated by 
a third region in which a non zero solution is instead 
stable. However the situation is even more complicated 
since, unlike the case Jo > 0, the non monotonicity of 
(3J^^^Xo{pJo,PJ^^^rn^^'>) refiects also in the fact that 
the self-consistent Eq. for E =F may have more 

solutions of the kind ±m^^\ ±m' {F) , . . . which are still 
stable with respect to the stability condition ([M]), for 
h = 0. We face in fact here the problem to compare more 
stable solutions. According to Eq. ([^5)1 . in the thermo- 
dynamic limit, among all the possible stable solutions, 
only m*^^\ the solution that minimizes L^^\ survives, 
whereas the not leading ones play the role of metastable 
states. This kind of scenario, which includes also finite 
jump discontinuities, has been besides observed in the 
context of small-world neural networks in plj where we 
even observe some analogy in the used formalism, at least 
for the simplest case of one binary pattern. 



8 



From Eqs. ([38]) and ([39)1 we see that the signs of 
the Landau coefficients a^^\ c^^\ . . ., are func- 

tions of /3 and Jq only. Given Jq < 0, the most im- 
portant quantity that features the non monotonicity of 

minimum value of /3 over 

which b^^^ becomes negative: 



6^^^ < 0, p>(3: 



(F) 



(58) 



The equation for as a function of Jq, defines a 

point where 6'^^-' = 0. If Jo < 0, the most general equa- 
tion for a generic critical temperature is no longer given 
by Eq. ((44|) . In fact, in general, a critical temperature 
now is any temperature where the stable and leading so- 
lution m^^^ may have a singular behavior, also with finite 
jumps between two non zero values. 

There are some simplification when for the Landau co- 
efficient c*^^\ we have c*-^' > 0, or at least c^^' > out of 
the P region In this situation in fact, from Eq. (|36p we 
see that a(^', b'-^^ and c^^) are the only relevant terms 
for the critical behavior of the system and - for small val- 
ues of m^^-* - we can again apply the Landau theory, this 
time for the so called model [20]. In such a case, for 
the solution m^^^ we have 



m 



(F) 



[^{b^'^y-^ai^^ci^)-bi^)y if, 



< 0, 



or 



a(F)>o and ^(f) < -4\/ (59) 



whereas 

m(J') = 0, if. 



/ (F) (F) 

a(F) > and b^^^ > -4^^-^. (60) 

From Eqs. and we see that, if 6^^) > 0, we 

have a second-order phase transition and Eqs. (|44)l - ([57|) 
are recovered with Eq. ([S^]) becoming the second of Eqs. 
PS)) for small and negative values of a^^^. However, from 
Eq. ([55)) we see that, if 6*^^^ is sufficiently negative, we 
have a first-order phase transition which, for small values 
of a^^\ gives 



6(F) 



1 - 



a(F)c(F) 



, if, 



a(F) < and 6^^) < g, or 



a(F) > and 6(f) < _4\/^-^. (61) 



From Eq. ^ we see that the line 6(f) = -4a/o(f)c(f)/3 
with a(F) > establishes a line of first-order transitions 
through which m(F) changes discontinuously from zero 
to 



Am(F) = 



3a(F)\^ 

c(F) 



(62) 



The point a(F) = 6(f) = is a tricritical point where 
the second and first-order transition lines meet. If we 
approach the tricritical point along the line 6(f) = we 
get the critical indices a — 1/2, a' = 0, /3 = 1/4, 7 = 
y — 1 and (5 = 5. However, this critical behavior along 
the line 6(f) = has not a great practical interest since 
from Eq. ([55[) we see that it is not possible to keep 6(f) 
constant and zero as the temperature varies. Finally, 
we point out that, even if c(f) > 0, when the transition 
is of the first-order, Eqs. ([59[) and ([6T|) hold only for 
6(F), and then a(F), sufficiently small, since only in such 
a case the finite discontinuity of 771(F) is small and then 
the truncation of the Landau free energy term ■0(f) to a 
finite order meaningful. Note that this question implies 
also that we cannot establish a simple and general rule to 
determine the critical temperature of a first-order phase 
transition (we will return soon on this point). 

When c(F) < 0, the Landau theory of the model 
cannot be of course applied. However, as in the case 
Jq > 0, even if the sign of c(F) cannot be a pri- 
ori established, for the convexity of the function fo 
with respect to f3h, the sum of the six-th term with 

A(/3/o) (/3J(^^\/3j(F)m(F)^, in Eq. §^ must go nec- 
essarily to +00 for m(F) — > 00 and a qualitative similar 
behavior of the mP model is expected. In general, when 
Jo < 0, the exact results are limited to the following ones. 

From now on, if not otherwise explicitly said, we shall 
reserve the name critical temperature, whose inverse 
value of P we still indicate with (3c , to any tempera- 
ture on the boundary of a P region (through which to(f) 
passes from to a non zero value, continuously or not). 
For each critical temperature, depending on the value of 
/3; , we have three possible scenario of phase transitions: 



4= first order, 

<^ tricritical point, 

^ second order. 



(63) 



Note that, according to our definition of critical tem- 
perature, the critical behavior described by Eqs. ([59[) - 
([62[) represents a particular case of the general scenario 
expressed by Eqs. ([63]) . We see also that, in general, 
when 6(F) < 0, approaching the tricritical point, for the 
critical exponent (3 we have (3 < 1/4. 

In the case in which pi^^ corresponds to a second-order 
phase transition, or in the case in which a(F) < out of 
the P region (at least immediately near the critical tem- 
perature), pi^^ can be exactly calculated by Eq. ([44[) . 
When we are not in such cases, the only exact way to 
determine the critical temperature is to find the full so- 
lution for m(F) which consists in looking numerically for 
all the possible solutions of Eq. ([T5[) and - among those 
satisfying the stability condition ([34l) - selecting the one 
that gives the minimum value of L(f) . 
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C. Level of accuracy of the method 

In the P region, Eqs. (|15ll33p are exact, whereas in the 
other regions provide an effective approximation whose 
level of accuracy depends on the details of the model. In 
particular, in the absence of frustration the method be- 
comes exact at any temperature in two important limits: 
in the limit c 0+, in the case of second-order phase 
transitions, due to a simple continuity argument; and in 
the limit c — *■ oo, due to the fact that in this case the 
system becomes a suitable fully connected model exactly 
described by the self-consistent equations (fTS]) (of course, 
when c — > oo, to have a finite critical temperature one 
has to renormalize the average of the coupling by c) . 

However, for any c > 0, off of the P region and in- 
finitely near the critical temperature, Eqs. (|15ll20p are 
able to give the exact critical behavior in the sense of 
the critical indices and, in the limit of low temperatures, 
Eqs. IjlSmop provide the exact percolation threshold. In 
general, as for the SK model, which can be seen as a 
particular model with Jq = 0, the level of accuracy is 
better for the F phase rather than for the SG one and 
this is particularly true for the free energy density f^'^\ 
Eq. pT|) . In fact, though the derivatives of /'^^ are 
expected to give a good qualitative and partly also a 
quantitative description of the system, /^^*^' itself can 
give wrong results when the SG phase at low tempera- 
tures is considered. We warn the reader that in a model 
with Jo = 0, and a symmetrical distribution dfi{Jij) with 
variance J, the method gives a ground state energy per 
site u^^'^\ which grows with c as m^^*^) ~ — Jc, whereas 
the correct result is expected to be u^^^^ ~ —J\fc p^ . 
As a consequence, in the SK model, in the limit /3 ^ cxd, 
the method gives a completely wrong result with an infi- 
nite energy. We stress however that the order parameters 
rrS^^ and vrS^'^^ , and then also the correlation functions, 
by construction, are exact in the zero temperature limit. 



derive exactly all the phase diagram. However, in our 
method, the solution with label F or SG are exact only 
in their own P region, i.e., the region where ra^'> = 
or m^^*^-* = 0, respectively. Unfortunately, according to 
what we have seen in Sec. IIIC, whereas the solution 
with label F is still a good approximation also out of the 
P region, in frustrated model (where the variance of d^i is 
large if compared to its first moment) the free energy of 
the solution with label SG becomes completely wrong at 
low temperatures. Therefore, we are not able to give in 
general the exact boundary between the solution with la- 
bel F and the solution with label SG, and in particular we 
are not able to give the physical frontier F/SG. However, 
within some limitations which we now prescribe, we are 
able to give the exact critical surface, i.e., the boundary 
with the P phase, establishing which one - in the ther- 
modynamic limit - of the two critical boundaries, P-F 
or P-SG, is stable (we will use here the more common 
expression "stable" instead of the expression "leading" ) , 
and to localize some regions of the phase diagram for 
which we can say exactly whether the stable solution is 
P, F, or SG. We will prove the stability of these solutions 
in Sec. IX. When for a region we are not able to discrim- 
inate between the solution with label F and the solution 
with label SG and they are both out of their own P re- 
gion, we will indicate such a region with the symbol "SG 
and/or F" (stressing in this way that in this region there 
may be also mixed phases and re-entrance phenomena). 

In Sec. IX we prove that there are four possible kind of 
phase diagrams that may occur according to the cases 
(1) (Jo > 0; do < 2, or Jo = oo), (2) (Jo > 0; 2 < 
do < oo), (3) (Jo < 0; Jo < 2, or do = oo), and (4) 
(Jo < 0; 2 < Jo < oo). The four kind of possible phase 
diagrams are schematically depicted in Figs. [l]|3]in the 
plane (T, X). 

1. Jo > 



D. Phase diagram 

The physical inverse critical temperature Pc of the ran- 
dom model is in general a non single- value function of X: 
Pc = Pc{X), where X represents symbolically the pa- 
rameters of the probability d/i for the couplings Jij, and 
the parameter c, the average connectivity (which is also a 
parameter of the probability distribution of the short-cut 
bonds) . The parameters of d/i can be expressed through 
the moments of dfi, and as they vary the probability dfi 
changes. For example, if d/i is a Gaussian distribution, 
as in the SK model, there are only two parameters given 
by the first and second moment. A concrete example for 
the one dimensional small-world model will be shown in 
Figs, (mi) and (UHl). 

In the thermodynamic limit, only one of the two so- 
lutions with label F or SG survives, and it is the so- 
lution having minimum free energy. In principle, were 
our method exact at all temperatures, we were able to 



As we have seen in Sec. IIIBl, if Jo > 0, for both the 
solution with label F and SG, we have one - and only one - 
critical temperature. In the following, to avoid confusion, 
it should be kept in mind the distinction between the 
physical Pc = Pc{X) and with E =F or SG. Pc 

satisfies the following rules. 

Case (1): If Jo < 2 and Jo is a finite range coupling, 
or else do = oo at least in a broad sense (see [lal), Pc{X) 
is a single- value function of X, and we have 

/3e = min{/3f),/3f«)}. (64) 

A schematic representation of this case is given in Fig. 1. 
Case (2): If instead 2 < Jo < oo we have 

pP>P.>pi'''\ if /3f^U/3f^ ^^'^ 

Notice in particular that the second line of Eq. (p5)) 
does not exclude that Pc{X) might be a non single- value 
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function of X. A schematic representation of this case is 
given in Fig. 2. 



2. Jo < 

As we have seen in Sec. IIIB2, if Jo < 0, for a suffi- 
ciently large connectivity c, the solution with label F has 
at least two separated P regions corresponding to two 
critical temperatures. Here we assume that the underly- 
ing lattice £o has only loops of even length so that, for 
example, triangular lattices are here excluded. Let us 
suppose to have for the solution with label F only two 
critical temperatures (the minimum number, if Jq < 0), 
and let be 



(66) 



In general we have the following scenario. 

Case (3): If do < 2 and Jo is a finite range coupling, 
or do = oo in a broad sense (see [l^), (ici{^^ is a single- 
value function of X. and satisfies Eq. ([51)) . The other 
critical inverse temperature /3ci(-X') is instead: either a 
two- value function of X. and we have 



/3c 



if 



in the cases (2) and (4) the SG phase is always limited 
below (in terms of temperatures) by the unstable P-SG 
surface coming from the solution with label SG (indicated 
as P-SG unst in Figs. 2 and 4). Finally, we stress that 
- under the hypothesis that Cq has only loops of even 
length - the stable P regions correspond always to the 
solution with label F. 

For 2 < do < oo, from the second line of Eqs. and 
(IM)) and from Eq. ([70)1 . we see that the method is not 
able to give the complete information about the P-SG 
boundary since we have only inequalities, not equalities. 
Furthermore, in these regions of the phase diagram the 
physical critical temperature in general may be a non 
single- value function of X. On the other hand, we have 
the important information that in these equations the 
inequalities between the physical Pc and /3c^^^ (the solu- 
tion with label SG) are always strict. As a consequence, 
we see that, when 2 < do < oo, in these regions the SG 
"magnetization" m^^^^ will always have a finite jump dis- 
continuity in crossing the surface given by (3c- In other 
words, along such a branch of the critical surface corre- 
sponding to the second line of Eqs. ([65]) and ([69)1 and 
Eq. (1701) ; have a first-order phase transition, indepen- 
dently of the fact that the phase transition corresponding 
to the /3c^'^'' surface is second-order, and independently 



(67) of the sign of Jo 



$ if /3^f^>/?f«\ 



(68) 



where $ in Eq. ([68]) means that if P^i > there 
is no stable boundary with the P region. A schematic 
representation of this case is given in Fig. 3. 

Case (4): If 2 < do < oo, /3c2 satisfies Eq. (105)) : 
whereas for p^i we have either 



cl 



(SG» 



if 



(69) 



if 3/^,1 ^/3,i>/?fG), if /3(p>/?fG), (70) 

where in Eq. (|69p we have introduced the symbol SG> to 
indicate that in general the stable P-SG surface is above 
(or below in terms of temperatures) the surface coming 

from the solution with label SG: /jf > Notice 
that, similarly to the case (3), we cannot exclude that 
Pel in Eq. (|70p be a non single-value function of X, as 



well as (3, 



(SG» 



in Eq. (|69p. A schematic representation 



of this case is given in Fig. 4. 

If more than two critical temperatures are present, the 
above scheme generalizes straightforwardly. 

Keeping our definition for the introduced symbol "SG 
and/or F" , we stress that: in all the fours cases the phases 
F and "SG and/or F" are exactly localized; in the cases 
(1) and (3) the phases P and SG are exactly localized; 



E. Generalizations 

The generalization to the cases in which the unper- 
turbed model has an Hamiltonian Hq involving couplings 
depending on the bond 6 S Fo is straightforward. In this 
case we have just to substitute everywhere in the for- 
mulae ([11])- dH]), J^^' with the set {Jq^^}. However, the 
critical behavior will be in general different and more 
complicated than that depicted in the Subsections IIIBl 
and IIIB2. In particular, even in the case in which all 
the couplings Jq^ are positive, we cannot assume that 
the Landau coefficient fc^^-* be positive so that, even in 
such a case, first-order phase transitions are in principle 
possible, as has been seen via Monte Carlo simulations 
in directed small- world models [23 |. 

As anticipated, our method can be generalized also to 
study possible antiferromagnetic phase transitions in the 
random model. There can be two kind of sources of anti- 
ferromagnetism: one due to a negative coupling Jo in the 
unperturbed model, the other due to random shortcuts 
Jij having a measure d/i with a negative average. 

In the first case, if for example the sublattice Cq is 
bipartite into two sublattices C^^^ and ^o'''', the unper- 
turbed model will have an antiferromagnetism described 
by two fields TOq"' and ml^\ Correspondingly, in the 
random model we will have to analyze two effective fields 
rn'"^ and to'''-' which will satisfy a set of two coupled self- 
consistent equations similar to Eqs. (|15p and involving 

the knowledge of m'j^^ and TOq^"* . More in general, we can 
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FIG. 1: (Color online) Phase diagram for the case (1): Jo > FIG. 3: (Color online) Phase diagram for the case (3): Jo < 

and do < 2 or do = oo in a broad sense. Here T is the and do < 2 or do = oo in a broad sense. T and X as in Fig. 

temperature while X represents symbolically the connectivity [T] 
c and the parameters of the probability distribution djj,. 





FIG. 2: (Color online) Phase diagram for the case (2): Jo > 
and 2 < do < oo. T and X as in Fig. [T] 



FIG. 4: (Color online) Phase diagram for the case (4): Jo < 
and 2 < do < cxD. T and X as in Fig. [l] 
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introduce the site-dependent solution moi to find corre- 
spondingly in a set of coupled equations (at most N), the 
eff'ective fields of the random model. 

In the second case, following [l^ we consider a lat- 
tice Cq which is composed of, say, p sublattices C'^\ 
V = 1, . . . ,p. Then, we build up the random model with 
the rule that any shortcut may connect only sites belong- 
ingto two different sublattices. Hence, as already done 
in [l3| for the generalized SK model, we introduce p effec- 
tive fields m^^^ which satisfy a system of p self-consistent 
equations involving the p fields ra"^^ and calculated in 
the p external fields J'^^^m*^") (note that here the symbol 
F stresses only the fact that the effective coupling must 
be calculated through Eq. (|16p ). 



IV. SMALL WORLD IN do = DIMENSION 

A. The Viana-Bray model 

As an immediate example, let us consider the Viana- 
Bray model It can be seen as the simplest small- 
world model in which N spins with no short-range cou- 
plings (here Jq — 0) are randomly connected by long- 
range connections J (possibly also random). Note that 
formally here £0 has dimension — Q. Since Jq = 0, for 
the unperturbed model we have 

-/3/o(0,/3/j) = log [2cosh(/3/i)], 



mo{0,Ph) = tanh(/3/i). 



Xo{PJo,l3h) = l-tanh2(/3/i)|/3^^o - 1, 



(71) 



It is interesting to check that the first and second deriva- 
tives oi xo in h = are null and negative, respectively. 
In fact we have 



and 



-|-Xo(0,/3/i) = -2tanh2(/3/i) 
pn 

X [l^tanh'ifSh)] \ph=o = 0, 



,Xo(0,/3M = -2[l-tanh2(/3;i)]' 



(72) 



4tanh2(/3/i) [l - tsmh^ {(3h)] \f3h=o = -2. (73) 

Applying these results to Eqs. (|15m9p we get immedi- 
ately the self-consistent equations for the F and the SG 
magnetizations 



TO^^-* = tanh 
TO^^^^ = tanh 



m'-^^c / d^tanh(/3J) 



(74) 



, (75) 



and the Viana-Bray critical surface 

d^tanh(/3f ) J) = 1, 
d^tanh2(/3f°)j) = 1. 



(76) 
(77) 



On choosing for dfi a measure having average and vari- 
ance scaling as 0(l/c), for c oc N, we recover the equa- 
tions for the Sherrington-Kirkpatrick model [1^ already 
derived in this form in [13 and dl. In these papers Eqs. 
(|74j|77p were derived by mapping the Viana-Bray model 
and, similarly, the Sherrington-Kirkpatrick model to the 
non random fully connected Ising model. In this sense 
it should be also clear that, at least for (3 < (3c and zero 
external field, in the thermodynamic limit, the connected 
correlation functions (of order k greater than 1) in the 
Sherrington-Kirkpatrick and in the Viana-Bray model 
are exactly zero. In fact, in the thermodynamic limit, 
the non random fully connected model can be exactly re- 
duced to a model of non interacting spins immersed in 
an effective medium so that among any two spins there 
is no correlation. Such a result is due to the fact that, in 
these models, all the N spins interact through the same 
coupling J /N , no matter how far apart they are, and the 
net effect of this is that in the thermodynamic limit the 
system becomes equivalent to a collection of N non in- 
teracting spins seeing only an effective external field (the 
medium) like in Eq. (j7ip with I3h replaced by (3Jm. 

For the measure pUj) our approximated Eq. ([Ti]) can 
be compared with the exact known equation that can 
be derived by using the Bethe-Peierls or the replica ap- 
proach and is given by (see for example [27j and refer- 
ences therein) 



,(F) 



9=0 



tanh(/3^ \[ «'(i/„)di/„. 



where the effective field H is determined by the following 
integral equation 



9=0 



(9-1)! 



5{H -Tt&nhr 



tanh (3 J x 



q-l q-1 

tanh(/3^ i/,„)]) Yl ■^{H„,)dH„ 

rn — 1 m— 1 



In the limit (3 ^ 00, Eqs. ^ and ([75]) give the fol- 
lowing size (normalized to 1) of the giant connected com- 
ponent 



,(SG) _ 



tanh(rn''^-'c), 
tanh(m^^^)c). 



(78) 
(79) 



These equations are not exact, however they succeed 
in giving the exact percolation threshold c = 1. In fact, 
concerning the equation ([75]) for the F phase, the exact 
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equation for m*^^) is (see for example [13] and references 
therein) 

l-m(F) = e-=™'", (80) 
which, in terms of the fmiction tanh, becomes 
2m(F) + (m(F))2 



2-m(F) + (to(f))2 



= tanh(m*-^-'c) 



so that Eqs. ([78)1 and ((80)) are equivalent at the order 
Oim^^'>). We see also that, as stated in the Sec. IIIC, 
Eqs. ([75]) and (|5D1) become equal in the limits c ^ and 
c oo. 



B. Gas of Dimers 



and - at least for Jq > - the equation for the critical 
temperature 



-/3f)j(^) =1. 



cosh(/?rvo 

As it will be clear soon, this model lies between the 
Viana-Bray model and the more complex = 1 dimen- 
sional chain small-world model, which will be analyzed 
in detail in the next section. Our major interest in this 
simpler gas of dimers small- world model is related to the 
fact that, in spite of its simplicity and do=0 dimension- 
ality - since the second derivative of xo may be positive 
when Jq is negative - according to the general result of 
Sec. IIIB, it is able to give rise to also multiple first- and 
second-order phase transitions. 



Let us consider for Cq a set of 2N spins coupled 
through a coupling Jq two by two. The expression "gas of 
dimers" stresses the fact that the dimers, i.e. the couples 
of spins, do not interact each other. As a consequence, 
the free energy, the magnetization, and the susceptibility 
of the unperturbed model can be immediately calculated. 
We have 

- PfoiPJo, f3h) = I log [2e^'» cosh{2(3h) + 2e-^''] , 



mo (/3 Jo, /3 ft.) 



g/JJo sinh(2/3ft) 



el^Jo cosh(2/3ft) -I- e~PJo ' 



2e'3"'o + 2 cosh(2/3ft) 
.0Jo cosh(2/3ft) + 6-/5-^0 



;2 1/9/1=0 



(81) 



cosh(/3Jo) ' 

Let us calculate also the second derivative of xo- From 

-|-Xo(/3Jo,/3ft) =4sinh(/3ft) 

-/3Jo _ 2e3/3Jo _ g/3Jo cosh(2/3ft) , , 

X -4-^, (82 

[el^Jo cosh(2/3/i) + e-'^-'o]^ 

sinh(/3Jo) + e^'^'^" 



we get 
52 



■Xo(/3Jo,/3/i)|/3h=o = -2- 



[cosh(/3Jo)]' 



We note that, as expected, the second derivative of xo in 
h — 0, for Jo > is always negative, whereas, for Jq < 
it becomes positive as soon as /?| Jo| > log(-\/2). 

By using the above equations, from Sec. Ill we get 
immediately the following self-consistent equation for the 
magnetizations 



tanh(2/3 J(^)m(^) + 2/3ft) 
1 -I- e-2/3Jo sech(2/3J(s)TO(s) + 2/3/i) ' 



V. SMALL WORLD IN do=l DIMENSION 

In this section we will analyze the case in which Cq is 
the c?o=l-dimensional chain with periodic boundary con- 
ditions (p.b.c). The corresponding small-world model 
with Hamiltonian in zero field has already been an- 
alyzed in by using the replica method. Here we will 
recover the results found in [5|] for f3c and will provide 
the self-consistent equations for the magnetizations m^^' 
and m^^^-' whose solution, as expected, turns out to be in 
good agreement with the corresponding solutions found 
in for c small and large (the latter when the frustra- 
tion is relatively small). It will be however rather evident 
how much the two methods differ in terms of simplicity 
and intuitive meaning. Furthermore, we will derive also 
an explicit expression for the two-points connected corre- 
lation function which, to the best of our knowledge, had 
not been published yet. Finally, we will analyze in the 
detail the completely novel scenario for the case Jo < 
which, as mentioned, produces multiple first- and second- 
order phase transitions. 

In order to apply the method of Sec. Ill we have to solve 
the one dimensional Ising model with p.b.c. immersed in 
an external field. The solution of this non random model 
is easy and well known (see for example [1^). If we 
indicate with Ai and A2 the two eigenvalues coming from 
the transfer matrix method, one has 



Ai,2 = e'^"'" cosh(^/i) ± [e^^^" 



from which it follows that, for the free energy density, the 
magnetization and the two-points connected correlation 
function, we have 

-/?/o(/3Jo,/3M = log(Ai), 



mo(/3Jo,/3/i) 



e'^-^" sinh(/3/i) 



(83) 
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Co(/3Jo,/?/»;||j-j|lo) 



def 



- [cr) 

2 



sm^{2tp) 



^XIK-Jllo 
Ai 



(84) 



where the phase ip is defined by 



cot(2v3) = e^^-'° sinh(/3/i), < ip < -, (85) 

and ||i — illg is the (euchdean) distance between i and j. 

Let us calculate xo and its first and second derivatives. 
From Eq. ((83|) we have 



g-/3Jo cosh(/3/i) 



[e2/5Josinh2(/3/i) +e-2/3Jo] 



(86) 



d(3h' 

[e-3/3Jo _ 2e'3"'» cosh2(/3/i) - e'^-'o] 
re2/3Josinh^(/3/i) + e-2/3-^<']* 



(87) 



.iXoiPJoiPh) = cosh(/3/i) 



a(/3/i)2 

^g-3/3.7o „ 2e/3"'o cosh2(/3/i) - e'^-'"] 
[e2/3Josinh2(/3/i) + e-Wo]* 



0(/3/i)'(88) 



From Eq. ([55]) we see that for Jq > and any /3 we 
have, for sufficiently small h, 

^p-y2Xo(/3Jo,/3/i) < cosh(/3/i) 



[1 - ief'M 



[e2/3Josinh2(/3/i) +e-2/3'A] 

whereas for Jq < we have 
92 



a(/3/i)- 



rXo(/3Jo,/3/i) > for e 



-4/3Jo 



>3. (90) 



We see therefore that, according to Sec. IIIB, when 
Jo < for (3\Jo\ > log(3)/4 = 0.1193... the Landau co- 
efficient b*-^-* is negative and we may have a first-order 
phase transition. 



From Eqs. ((T5)) and ((83)) . for the magnetizations m^^^ 
and m*^^*^-* at zero external field we have 



TO 



e'5^o'^'sinh(/3j(^)TO(^)) 



e^/'^''' sinh2(/3J(S)m(5^)) + e-^/J^' 



(91) 



From Eqs. (|44p and ([5^ we see that a solution to*^^^ 
becomes unstable at the inverse temperature (3c given 
by 

For Jo > the above equation gives the exact P-F and 
P-SG critical temperatures in agreement with Q. When 
Jo < - unless the transition be of second-order - Eq. 
([92| for E =F does not signal a phase transition. In 
general, as Jo < the P-F critical temperature must be 
determined by looking at all the stable solutions to*^^^ of 
the self-consistent equation ([9T|) and by choosing the one 
minimizing the effective free energy L^^\m) of Eq. (j22|) . 

Finally, for the two-point connected correlation func- 
tion, from Eqs. ([20]), ^ and dM]), we have 



C(^)(||i-j||J = sin2(2^(^))e-l'^-^Ho/«' 



where 



2^(^) = cot-i [e^^^^"' sinh(/3j(^)TO(^))' 



and the correlation length S}'^^ is given by perform- 
ing the effective substitutions /?Jo PJg and l3h — > 
/3j(^)m(^) in log(Ai/A2). 

Note that Co(/3Jo, /?/i) is even in (3h, so that C{—m) = 
C{m). Near the critical temperature we have 



sin(2v3(^^) = 1 - 



+ 0(to(^))4 



and 



log[tanh(/3J^"0] 



(S),^ (/3J(^)to(^))2 



4sinh(/3J(^ 



' -I- e3/34"') tanh(/3jf )) + e^'^^o"' _ e/^^"'] + 



According to the general result, Eqs. ([55)1 - (|57|) . we see that the correlation length remains finite at all tempera- 



15 



One second-order phase transition 



No phase transition 




(F) 

m 

(F) (F) 

Stab, or m and m =0 



FIG. 5: (Color online) Magnetization (thick solid line), and 
curves of stability (dashed and dot-dashed lines) for the case 
c = 0.5, Jo = 1 and J = 3/5/. Here = 1.687. 

One second-order phase transition 
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/ \ 
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_ (F) 
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stab, of m'"'^* 
stab, of m^'^'=0 
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1,1 



FIG. 6: (Color online) Magnetization (thick solid line), and 
curves of stability (dashed and dot-dashed lines) for the case 
c = 10, Jo = 0.25 and J = 1/c. Here Tc = 1.419. 



FIG. 7: (Color online) Magnetization (thick solid line), and 
curves of stability (dashed and dot-dashed lines) for the case 
c = 5, Jo = —1 and J = 1. Note that here m'^' = and the 
two curves of stability coincide everywhere. 

No phase transition 



0.4 - / 



(F) 

m 

stab, of m*'^' and m*^'=0 



3 4 5 6 7 

T 



FIG. 8: (Color online) Magnetization (thick solid line), and 
curves of stability (dashed and dot-dashed lines) for the case 
c = 5.828, Jo = -1.4 and J = 1. Note that here m'^^ = 
and the two curves of stability coincide everywhere. 



tures. 

In Figs. [51IT^ we plot the stable and leading magneti- 
zation m(F) (thick sohd line), xq [i3^^U^^\q^(3^^U^^'^ 

(dashed hne), and xo \ /3J^^^"i(f))^(f) J(f) 

(dot-dashed Une) for several cases obtained by solving 
Eq. (|9ip numerically with E=F. The stable and leading 
solution (the only one drawn) corresponds to the solu- 
tion that minimizes L(^)(to) (seeEq. In aU these 
examples we have chosen the measure (|10l) . Figs. [S]and 
[S] concern two cases with Jq > so that one and only one 
second-order phase transition is present. The input data 
of these two cases are the same as those analyzed numer- 
ically in Pi (note that in the model considered in Q , the 
long range coupling J is divided by c). As already stated 
in Sec. Ill, the self-consistent equations become exact in 
the limit c — > 0, for second-order phase transitions, and 
in the limit c — s- cx). Therefore, for the magnetization, 
by comparison with [5|, in Figs. [n]and[Sl where c is rela- 
tively small and big, respectively, we see good agreement 



also below the critical temperature. 

Figs. [TlfHl concern eight cases with Jo < 0. In these 
figures we plot also the line y = 1 , to make evident when 
the stability conditions for the solutions m^^^ — and 

to(I') ^ 0, which are given by xo [p^^'^ jfK^^P^'^U^'^^ < 

1 (dashed line), and xo (^P^^^ 4P , PJ^^^m^^^^P^^U^^') < 

1 (dot-dashed line), are satisfied, respectively. As ex- 
plained above, the critical behavior and the localization 
of the critical temperatures is more complicated when 
Jo < 0. In particular, given |Jo|, if c is not sufficiently 
high the solution m'^^^ — remains stable at all tem- 
peratures and if it is also a leading solution, no phase 
transition occurs. Let us consider Eq. ([5^ . For Jo < 
the Ihs of this equation has some maximum at a finite 
value f3 given by 
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T,,=3.641, two first-order phase transitions T^=1.820, first- and second-order phase transition 




4 5 

T 



FIG. 9: (Color online) Magnetization (thick solid line), and 
curves of stability (dashed and dot-dashed lines) for the case 
c = 5.5, Jo = -1, and J = 1. Here Td = 1.02 and Tc2 = 2.27. 



T,,=3.277, two first-order phase transitions 



1 1 1 1 1 1 1 1 1 




_ (F) 

m 

stab, of m*^' 
' stab, of m'^^'=0 




/ 

- / 

/ 


,1 


,1 


1,1,1, 





FIG. 10; (Color online) Magnetization (thick solid line), and 
curves of stability (dashed and dot-dashed lines) for the case 
c = 5.5, Jo = -0.9 and J = 1. Here Td = 0.85 and Tc2 = 
2.78. 

def 

where r — | Jo|/J, and we have introduced 

S{r) \/l + r"^ — r. 

Hence, we see that a sufficient condition for the solution 
r7T,(^) = to become unstable is that be 

Note that the above represents only a condition for the 
instability of the solution m^^^ = 0, but the true solution 
is the one that is both stable and leading. In fact, when 
Jo < 0, a phase transition in general may be present also 
when Eq. ([93|l is not satisfied and, correspondingly the 
possible critical temperatures will be not determined by 
Eq. (IMl). 

In Fig. [7]we report a case with J = 1, Jg = — 1 (r = 1) 
and a relatively low value of c, c = 5, so that no phase 
transition is present. Similarly, in Fig. [8]we report again 
a case in which no phase transition is present due to the 



FIG. 11: (Color online) Magnetization (thick solid line), and 
curves of stability (dashed and dot-dashed lines) for the case 
c = 6, Jo = -0.5 and J = 1. Here Td = 0.35 and Tc2 = 4.87. 
Tc2 corresponds to a second-order phase transition. 



T^=0.728, first- and second- order phase transitions 
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T 

FIG. 12; (Color online) Magnetization (thick solid line), and 
curves of stability (dashed and dot-dashed lines) for the case 
c = 4, Jo = -0.2 and J = 2. Here Td = 0.22 and Tc2 = 7.55. 
Tc2 corresponds to a second-order phase transition. 



fact that here r is relatively big, r 1.1. It is interesting 
to observe that for r = 1 Eq. requires a value of c 
greater than the limit value c = 3 + 2\/2 = 5.8284.... In 
both Figs. [QlandflOlwe report a case in which Eq. (|93p is 
still not satisfied, but nevertheless two first-order phase 
transitions are present. In both Figs. [TT] and [T^ we have 
one first- and one second-order phase transition. In both 
Figs. fT3l and [Ml we have two second-order phase transi- 
tions. As anticipated in Sec. IIIB, we note that in Figs. 
[9lfT2l i.e., the cases in which there is at least one ffi-st- 
order phase transition, there are always regions where 
both the solutions, m'^^ = and m^^^ ^ 0, are simulta- 
neously stable but only one solution is leading (the one 
drawn); whereas in Figs. fT3landfT4l as it was in Figs. [5] 
and [HI since we have only second-order phase transitions, 
the stability condition turns out to be a necessary and 
sufficient condition for determining the leading solution 
and the critical temperature can be determined also by 
Eq. (HH) with S =F. 



17 



T,=2.185, two second-order phase transitions Spin glass second-order phase transition 




FIG. 13: (Color online) Magnetization (thick solid line), and 
curves of stability (dashed and dot-dashed lines) for the case 
c = 1.6, Jo = -0.6 and J = 7. Here Tci = 2.58 and Tc2 = 
7.55. 



7,^=1.820. two second-order phase transitions 
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FIG. 15: (Color online) Spin glass order parameter (solid 
line), and curves of stability (dashed and dot-dashed lines) 
for the case c = 0.5, Jo = 1 and J = 3/5/c. Here Tc — 1.130. 

Spin glass second-order phase transition 



0.8 - 



0.6 - 



0.4 - 



0.2 - 




FIG. 14: (Color online) Magnetization (thick solid line), and 
curves of stability (dashed and dot-dashed lines) for the case 
c = 1.4, Jo = -0.5 and J = 10. Here Td = 3.00 and Tc2 = 
9.34. 



FIG. 16: (Color online) Spin glass order parameter (solid 
line), and curves of stability (dashed and dot-dashed lines) 
for the case c = 10, Jo = 0.25 and J = 1/c. Here Tc = 0.424. 



In the top of Figs. [9lfT4lwe write the discriminant tem- 
perature = 4|Jo|/log(3) below which a phase transi- 
tion (if any) may be first-order [see Eq. and Eqs. 

FinaUy in Figs. [15] and [16] we plot the 
spin glass order parameter to'^*^' (solid line), 
Xo (^^^°^Jo'^°\o)/5(SG)j(SG) (dashed line), and 

Xo (/?(SG)j(SG)^^j(SG)^(SG)^^(SG)j(SG) (dot-dashed 

line) obtained by solving Eq. (j9ip numerically with 
E=SG. In these two examples we have chosen the mea- 
sure (jll[) and, for c, \J\, and Jo, we have considered the 
same parameters of Figs. [5] and [6] of the ferromagnetic 
case. As anticipated in Sec. IIIB, due to the fact that 
the effective coupling Jg^*^' is positive, there is only a 
second-order phase transition, the stability condition 
turns out to be a necessary and sufficient condition 
for determining the leading solution, and the critical 
temperature can be determined also by Eq. (|92p with 



E =SG. 

Note that, unlike the P-F critical surface, the P-SG 
critical surface does not depend on the parameter p enter- 
ing in Eq. pTjl . For the reciprocal stability between the 
P-F and the P-SG critical surface we remind the reader 
to the general rules of Sec. HID (see cases (1) and (3)) 
which, for Jq > 0, reduce to the results reported in Sec. 
6.1 of the Ref. [B]. Here we stress just that, if Jq > 0, for 
p < 0.5, only the P-SG transition is possible. However, 
when Jo < and c is not sufficiently large, the SG phase 
may be the only stable phase even when p = 1. In fact, 
although when Jo < the solution m*-^' may have two 
P-F critical temperatures, in general, if the P-SG tem- 
perature is between these, we cannot exclude that the 
solution m^^*^^ starts to be the leading solution at suf- 
ficiently low temperatures. In Figs. [T7] and [TH] on the 
plane (T, c), we plot the phase diagrams corresponding 
to the cases of Figs. [TSIandfni respectively. These phase 
diagrams are obtained by solving Eq. supposing that 
here, as in the cases of Figs. [T3]and[T3] where c = 1.4 
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Phase diagram 
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FIG. 17: (Color online) Phase diagram for the case considered 
in Fig. [13] with the measure of Eq. (fTU)) . 



Phase diagram 



Tio 




P-F 
P-SG 
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FIG. 18: (Color online) Phase diagram for the case considered 
in Fig. [14] with the measure of Eq. (|10|) . 



and c = 0.5, respectively, the P-F transition is always 
second-order. We plan to investigate in more detail the 
phase diagram in future works. 



VI. SMALL- WORLD SPHERICAL MODEL IN 
ARBITRARY DIMENSION do 

In this section we will analyze the case in which 
the unperturbed model is the spherical model built up 
over a do-dimensional lattice Co (see [11] and references 
therein) (sj. In this case the cr's are continuous "spin" 
variables ranging in the interval (— oo,cxd) subjected to 
the sole constrain 7^,-^^ f^f — N, however our theorems 
and formalism can be applied as well and give results 
that, within the same limitations prescribed in Sec. Ill, 
are exact. 

Following for the unperturbed model we have 



1 



l3fo{pJo,Ph) = - log 



TT 



+ (l){pJo,l3h,z), 



moiPJo,Ph) 



(3h 
2/3 Joz' 



(94) 



where 

0(/3Jo,/3/i,z) 

9i^) 



PJodo + pJoz 



1 , , i(3hy 



(27r)'^o JO JO 



. duj, 



do 



X log [do + 2 - cos(tJi) - ... - cos(cJdJ] , 

and z — z{/3Jq, (3h) is the (unique) solution of the equa- 
tion dz<j) {(iJoiPh, z) — 0: 



/3Jo 



1 '(-\ 
2^ 



(95) 



(96) 



(97) 



Jq {it) being the usual Bessel function whose behavior 
for large t is given by 



4/3Joz2 

from which it follows the equation for toq 

The derivative g' can in turn be expressed as 
g'{z) - / [Jo {it)t' dt, 



Jo (it) 



(27ri)3 



1 




The critical behavior of the unperturbed system de- 
pends on the values of g'{z) and g"{z) near z — 0. It 
turns out that for do < 2 one has g'{0) = oo and there 
is no spontaneous magnetization, whereas for do > 2 one 
has g'{0) < oo and aX h = the unperturbed system 
undergoes a second-order phase transition with magneti- 
zation given by Eq. which, for /3 above /3cOi becomes 



TOo(/3Jo,0) 



where the inverse critical temperature (3co is given by 

PcoJo = la' (0) . 

Furthermore, it turns out that for do < 4 one has 
g"{0) = oo, whereas for do > 4 one has ^"(0) < oo. 
This reflects on the critical exponents a, 7 and d, which 
take the classical mean-field values only for do > 4. 

According to Sec. Ill, to solve the random model - 
for simplicity - at zero external field, we have to per- 
form the effective substitutions /SJq —* fiJ^ and I3h 

m^^-* in the above equations. From Eqs. (I94p . (|95p 
and (|96p . we get immediately: 



2/3jf) 
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the equations for inverse critical temperature P\ 



2^ \20hP^ 
and the magnetizations mP"^ 



0, /3 < /3f ^ > 0. 



Note that, as it must be from the general result of Sec. 
IIIB, unlike the unperturbed model, as soon as the con- 
nectivity c is not zero, Eq. (^5)) has always a finite so- 
lution /3c , independently on the dimension do- In fact, 
one has a finite temperature second-order phase transi- 
tion even for do 0+ where from Eq. (|97p we have 



[do = 0) 



so that the equations for the critical temperature 
become 



1, 



(do = 0) 



which, as expected, coincide with Eqs. (I76l) and (|77p of 
the Viana-Bray model. 

(98) Similarly, unlike the unperturbed model, in the ran- 
dom model all the critical exponents take the classical 
mean- field values, independently on the dimension dg. 
In the specific case of the spherical model, this behavior 
is due to the fact that g'{z) and g"{z) can be singular 
only at z = but, as soon as the connectivity c is not 

(99) zero, there is an effective external field /3J'^^^m*^^^ so that 



is not zero. For the critical behavior, the dependence 



on the dimension do reflects only in the coefficients, not 
on the critical exponents. In particular, concerning the 
argument of the square root of the rhs of Eq. by 
expanding in the reduced temperature t^^^ , for li*^^^ | <C 1 
we have 



where 



(S) t(S) 







(100) 



,(2/?fvr)^ 



1 - tanh^ 



(F) 



))/3f)j,,, -c J dA^(J.,, )tanh(/3f)j,j) 



^(SG) 



2c 



, tanh(/3f^)jo)/?f^)jo , J 



l + tanh2(/3f°^Jo) 



(2/3r"V^^"0 
d^l{AJ) (l - tanh'(/3(SG) tanh(/3f j^^^.)^(SG) 

tanh(/3^'^'^)jo)/3(^'^^Jo 



4c / d^^{J^J)tSinh^{Pi^'^^J^J) 



l + tanh2(/3f^)jo) /3f«V(SG) 



so that from Eqs. (l99|) and ()100|) for the critical behavior 
of the magnetizations we get explicitly the mean-field 
behavior: 



0, i(^) > 0. 



VII. MAPPING TO NON RANDOM MODELS 

In Sec. VIII we will derive the main result presented in 
Sec. III. To this aim in the next Sec. VIIA we will recall 
the general mapping between a random model built up 
over a given graph and a non random one built up over 
the same graph, whereas in Sec. VIIB we will generalize 
this mapping to random models built up over random 
graphs. We point out that the mapping does not consist 
in a sort of annealed approximation. 
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A. Random Models defined on Quenched Graphs 

Let us consider the following random model. Given a 
graph g, which can be determined through the adjacency 
matrix for shortness also indicate by g = {gt}, with = 
0, 1, 6 being a bond, let us indicate with Tg the set of the 
bonds b of g and let us define over Tg the Hamiltonian 

H {{ai};{Jb}) "^^^ - ^ Jb(Jit<^jt -^/ijtTi (101) 

beFg i 

where Jt is the random coupling at the bond 5, and 
tTijj , (7i^ are the Ising variables at the end-points of b. The 
free energy F and the physics are defined as in Sec. II 
by Eqs. ®-®: 

- pF 1/ J dV ({ JJ) log (Z ({J,})) , (102) 

W'= Jdri{Jb}){Oy, 1 = 1,2 (103) 

where dV {{Jt}) is a product measure over all the possible 
bonds b given in terms of normalized measures dfib > 
(we are considering a general measure dfit allowing also 
for a possible dependence on the bonds) 



dV{{.h}) 



Y[ dfib {Jb) ■ 



rf/ifc ( Jfc) = 1, (104) 



where Ffuu stands for the set of bonds of the fully con- 
nected graph. As in Sec. II, we will indicate a generic 
correlation function, connected or not, by C with under- 
stood indices ii,...,ik all different, see Eqs. and 

In the following, given an arbitrary vertex i of g, we 
will consider as first neighbors j of i only those vertices 
for which J d^i,j{Jij)Jij or J diJ,ij{Jij)Jfj are at least 
0{i/N). Note that we can always neglect couplings hav- 
ing lower averages. We will indicate with D{Tg) the av- 
erage number of first neighbors of the graph g. For a 
d-dimensional lattice, D{Tg) = 2d— 1, for a Bethe lattice 
of coordination number q, D(Tg) = q — 1, and for long 
range models, D{Tg) oc N. We will exploit in particular 
the fact that D{Tco U Tfun) = i:'(rfun) oc N. 

Given a random model defined trough Eqs. (|10im04p . 
we define, on the same set of bonds Tg, its related Ising 
model trough the following Ising Hamiltonian 

Hi{{<J^};{4'^}) =^ - j(^V,,aj,-^M,a05) 

&Gr„ i 



In the following a suffix / over quantities such as Hj, Fj, 
fi, gi, etc. . . , or /3i^\ etc. . . , will be referred to the 
related Ising system with Hamiltonian (jlOSp . 

We can always split the free energy of the random sys- 
tem with N spins as follows 

-(3F = Y f dub log [cosh(/?Jfc)] + 

b 

log [2 cosh(/3/ii)] + 4>, (108) 

i 

(j) being the high temperature part of the free energy. Let 
Lf be the density of cj) in the thermodynamic limit 



09 lim (b/N. 

N^oo 



(109) 



Let us indicate with ipi the high temperature part of 
the free energy density of the related Ising model defined 
through Eqs. (|105m07p . As is known, (^/ can be ex- 
pressed in terms of the quantities zi, = tanh(/9Jj^^^) and 
Zi = tanh{phi), i.e., the parameters of the high temper- 
ature expansion: 

ipi = ipi ({tanh(/3J,(^))}; {tanh(/3/i,)}) • (HO) 

The related Ising model is defined by a set of, typically 
few, independent couplings {jj,^^}, trough Eqs. (|106I 
I107P and, for /i^ = 0, i = 1, . . . ,N, its critical surface will 
be determined by the solutions of an equation, possibly 
vectorial, Gj ({tanh(/3jy^)}) = 0. 

In [Tsj we have proved the following mapping. 

Let /3c^^'' and /Sl^^^^"* be respectively solutions of the 
two equations 

G/ (^{y"dM6tanh2(/3(SG)j^)}^ =0, (111) 
Gi ({ I dM6tanh(/3f /AF)jfc)}) = 0. (112) 



Asymptotically, at sufficiently high dimensions D{Tg), 
the critical inverse temperature of the spin glass model 
/3c is given by 



/3. = min{/3f«),/3f/^^)}; 



(113) 



and in the paramagnetic phase for D{Tg) > 2 the follow- 
ing mapping holds 



(P - (Peff 



C~G, 



eff 



c 



= o 



( 1 



(114) 



where the Ising couplings J^^ have non random values 
such that y b,b' € Tg 

= jj,^^ if d^b' = d^ib. 



(106) 



feff = j^^""^ rfMfctanh'(/3Jb)} ) , (115) 



and 
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where 



2, if <fi{{JdfibUnh''{pJb)})> 
2\ipi ({/ dM6tanh(/3Jb)}) I, 
1, if (/?/({/ dAi6tanh^(/3Jb)}) < 
2|^/ ({/ dfibtanhiPJb)}) I, 



(117) 



and S^F/AF or SG, for /=1 or 2, respectively. 

In the hmit D(rg) oo and /i^ — 0, i = 1, . . . , iV, 
Eqs. pimiiyp . give the exact free energy and corre- 
lation functions in the paramagnetic phase (P); the ex- 
act critical paramagnetic-spin glass (P-SG), /3c , and 



(F/AF) 



surfaces, whose 



paramagnetic- F/AF (P-F/AF), j3, 
reciprocal stability depends on which of the two ones has 
higher temperature. In the case of a measure dfi not 
depending on the bond b, the suffix F and AF stand 
for ferromagnetic and antiferromagnetic, respectively. In 
the general case, such a distinction is possible only in 
the positive and negative sectors in the space of the pa- 
rameters of the probability distribution, {J dfibJb > 0} 
and {J d^bJb < 0}, respectively, whereas, for the other 
sectors, we use the symbol F/AF only to stress that the 
transition is not P-SG. 

It is not difficult to see that, when the measure d/i 
does not depend on the specific bond b, i.e., if d/if, = 
d^J.b' yb,b' £ Fg, in the P region Eqs. (I11HI117P lead to 
the following exact limit for tp and C ^ISj 

lim ip= lim C = 0, for (3 < (3c, (118) 

therefore, the basic role of Eqs. (|114mi7|) . is to show 
how, in the limit D{Tg) ~* oo, tp and C approach zero 
and which are their singularities. In particular this proves 
that for all the (random) infinite dimensional models and 
any disorder non bond-dependent, the critical exponent 
a' for the specific heat has the mean- field classical value, 
a' — 0, and that the correlation functions (with differ- 
ent indices) above the critical temperature are exactly 
zero. We point out however that, when the measure dfXb 
depends explicitly on the bond b, Eq. (jllSp in general 
does not hold [32| . In fact, when the measure dfib is 
bond-dependent, the symmetry expressed by Eq. (jllSp 
is broken since the bonds are no longer equivalent. As we 
will see in the next section, in small-world models with 
an underlying lattice £o having do < 2, even if Eq. (|118p 
may still holds for tp, the symmetry is broken for C since 
the direction(s) of the axis(es) of Cq is(are) now favored 
direction(s). Yet, if 2 < dp < oo, the symmetry (IllSp for 
if is broken as well. 

The analytic continuation of Eqs. pi4j|117p to /? > 
/3c and/or for h ^ provide certain estimations which 
are expected to be qualitatively good. In general such 
estimations are not exact, and this is particularly evident 
for the free energy density of the SG phase. However, 
the analytic continuation for the other quantities gives 
a good qualitative result and provide the exact critical 
behavior (in the sense of the critical indices) and the 
exact percolation threshold. 



For models defined over graphs satisfying a weak def- 
inition of infinite dimensionality, as happens on a Bethe 
lattice with coordination number g > 2, a more general 
mapping has been established 15|. In this case, all the 
above equations - along the critical surface (at least) - 
still hold exactly in the thermodynamic limit, where we 
can set effectively D{Tg) = oo. However, for the aims of 
this paper we do not need here to consider this general- 
ization of the mapping. 

We have yet to make an important comment about 
Eqs. (^5]) . and (^5]) . concerning the evaluation of 
a correlation function in the SG phase here for a ran- 
dom system with Jq = (for the moment being). In fact 
Eq. (|116p . for both a normal and a quadratic correla- 
tion function C*^^^ or C'^\ has a factor 1/2 not entering 
in the physical Eqs. ([23, (EH) and The difference 

is just due to an artefact of the map pin g that separates 
the Gibbs state into two pure states [18| not only in the 
F case, but also in the SG case. In fact, let us con- 
sider the correlation functions of order fc = 1, that is, 
C(i) = and C(2) = = (cri)2. We see that, for 
C(i), Eq. (inSl) in the SG phase gives Ct^' = m'-^'^'^ /2. 
On the other hand, for any non zero solution m'^^^ of 
the self-consistent Eq. ([15]), there exists another solu- 
tion —m^^'^\ and both the solutions have 1/2 probabil- 
ity to be realized in the random model. Since the SG 
phase is expected to be the phase characterized by having 
Qea 7^ and (ai) — 0, we see that if we introduce both 
the solutions m'^^^ and —m^^'^\ we get (cti) = in the 
SG phase. Similarly, for C^^), Eq. in the SG phase 



gives 



) /2, which at zero temperature gives 



1/2, whereas a completely frozen state with qea = 1 is 
expected. Again, we recover the expected physical qea 
by using both the solutions m*-^*^^ and —m^^'^\ Repeat- 
ing a similar argument for any correlation function of 
order fc, and recalling that for fc even (odd) the correla- 
tion function is an even (odd) function of the external 
magnetic field h, we arrive at Eqs. (pS)) and 



B. Random Models defined on Unconstrained 
Random Grapiis 

Let us consider now more general random models in 
which the source of the randomness comes from both the 
randomness of the couplings and the randomness of the 
graph. Given an ensemble of graphs g € Q distributed 
with some distribution P{g), let us define 

Hg {{<J^}■, {Jb}) =^ - Jb<J^i,<^j, -hYa^ 

beFg i 

= - XI 9bJb(yi,(Tj^-h'Yi[llQ) 

berf„ii « 

The free energy F and the physics are now given by 

- I3F J2 P^a) [ dV {{Jb}) log [Zg {{Jb})) , 
geg 
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and similarly for (O)', I = 1,2. Here Zg {{Ji,}) is the 
partition function of the quenched system onto the graph 
realization g with bonds in Tg 

and dV {{Jb}) is again a product measure over all the 
possible bonds b given as defined in Eq. p04p . Note 
that the bond-variables {gb} are independent from the 
coupling- variables {Jb}- 

For unconstrained random graphs, or for random 
graphs having a number of constrains that grows suffi- 
ciently slowly with N, the probability P{g), for large N, 
factorizes as 

pi9)= n 

In such a case we can exploit the mapping we have pre- 
viously seen for models over quenched graphs as follows. 
Let us define the effective coupling Jb: 

Jb = Jb ■ gb, 

correspondingly, since the random variables Jb and gb are 
independent we have 

dp.b{Jb) = dnb{Jb) -Pbigb), 
with the sum rule 

/ djib{Jb)f{Jb\gb)^ ^ Pb{gb) d^ibiJb)f{Jb]gb)- 

•' 36=0,1 

As a consequence, if we define the following global mea- 
sure 

dV [{Jb}) = Pig) ■ dV {{Jb}) = n dfibiJb), 

6ert„ii 

we see that the mapping of the previous section can be 
applied as we had a single effective graph Tp given by 

Tp {b e Tfuu : Pb{gb = 1) 7^ 0}, 
in fact we have 

-f3F^ JdV ({J,}) log (Zp {{Jb})) , 

where Zp is the partition function of the model with 
Hamiltonian Hp given by 

Hp{{a,};{Jb}) - Jb'^H<^j.~hY,'^,- (120) 



VIII. DERIVATION OF THE 
SELF-CONSISTENT EQUATIONS 

By using the above results, we are now able to derive 
easily Eqs. (|15ll23p . Sometimes to indicate a bond b we 
will use the symbol {i,j), or more shortly ij. 

It is convenient to look formally at the coupling Jq also 
as a random coupling with distribution 

d/.o(4)/d4 = 5( J'^ - Jo). (121) 
Let us rewrite explicitly the Hamiltonian ((2]) as follows 
iJc = - ^ {Jo + Ci-jJi-j) aia-j 

~ X/ CijJijaiCTj - h'Y^cFi, (122) 

and let us introduce the random variables J^, g'f, and J^, 
where 

, def ( Jo + CbJb, be To, 
'''' - \Jb, bi To, 



, def ( 1, be To, 
^ \cb, bi Fo, 

and 

•Jb — 'Jb ' 9b- 

Taking into account that the random variable Jg -I- Cij Jij , 
up to terms 0{1/N), is distributed according to dfio{Jo), 
the independent random variables and 5^ have distri- 
butions dfi'ij and respectively given by 

'^^'b^^^>-\d^l{J',), biVo, 

and 

"'^''^^[pia'b), hiu, 

where the measures dfi and p are those of the model 
introduced in Sec. II. As a consequence, Eq. (|122p can 
be cast in the form of Eq. (|120p with the measure 

^''^ ^' \ dfi{Jl,)p{g'f,), biVo. 

Finally, since Pb{gb) 7^ for any b e Ffuii, we have also 

Tp = Ff„ii, (124) 

and due to the fact that -D(Ffuii) cx N, in the thermody- 
namic limit the mapping becomes exact. 
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According to Eqs. (|105m07|) . the related Ising model 
of our small-world model has the following Hamiltonian 

with two free couplings: Jq'^\ for Fq, and J^^\ for Ffuu 



Hr 



-J, 



{I) ST 





(i,i)ero 



(125) 



After solving this Ising (/) model the mapping allows 
us to come back to the random model by performing 
simultaneously for any b G Tf^n the reverse substitutions 



tanh (^PJ^^' 



d/iUJ6)tanh'(/34), (126) 



where / = 1, 2 for E = F or SG solution, respectively. 
Since the couplings Jq^^ and J^^^ are arbitrary, we find it 
convenient to renormalize J^^'^ as J^^^ /N and at the end 
of the calculation to put again j(-f) instead of j(^) /N. 
Note that for the mapping nothing changes if we do not 
make this substitution; the choice to use j'^^'^ /N instead 
of j'^'' is merely due to a formal convenience, since in 
this way the calculations are presented in a more stan- 
dard and physically understandable form. In fact, ac- 
cording to Eqs. p23p and (|126p what matters after solv- 
ing the related Ising model with j'^^^ /N instead of J^^-* 
is that, once for S=F and once for I]=SG, we perform 
- simultaneously in the two couplings - the following re- 
verse mapping transformations {I — 1,2 for S = F or SG, 
respectively): 

tanh [pj^^'^/N^ j djl (jyj tanh' (/3Jy) (127) 

for (i, j) ^ Fo, and 

tanh (/3J^^') j djl (jy) tanh' (/3Jy) , (128) 

for e Fq. 

Explicitly, by applying Eqs. (fT^ . ^ and ([Hi]), the 
transformations (|127p and (|128|) become, respectively 



and 



13 jW ^ 13 ji^) 



(129) 



(130) 



where we have made use of the definitions (fTBl) - P^ in- 
troduced in Sec. III. 

Let us now solve the related Ising model. We have to 
evaluate the following partition function 

In the following we will suppose that J W (and then J(^)) 
is positive. The derivation for j'^-'^ (and then j'--^^) neg- 
ative differs from the other derivation just for a rotation 



of 7r/2 in the complex m-plane, and leads to the same 
result one can obtain by analytically continue the equa- 
tions derived for J^^' > to the region j'^' < 0. 

By using the Gaussian transformation wc can rewrite 
Zj as 



Zi ^ cn'Y e' 



where c^r is a normalization constant 



Cat 



2^ ' 



and, in the exponent of Eq. p3ip . we have again ne- 
glected terms of order 0(1). For finite N we can exchange 
the integral and the sum over the ct's. By using the def- 
inition of the unperturbed model with Hamiltonian iJg, 
Eq. Ill), whose free energy density, for given /3Jo and /3ft., 
is indicated with /o(/3 Jo, /3ft), we arrive at 



Zi = cn 



dm e-^^(™) 



(132) 



where we have introduced the function 

L(r7i) = ^j(^)m2H-^/o(/3J($^\/3j(-^)m-H/3ft) . (133) 

By using dph Ph{PJo,Ph) = -mo{/3Jo, ph), and 
di3h mo{pJo,Ph) = xo(/3Jo,/3ft) we get 



L' 



(ml \m — mo 



L" 



(m) = /3J(^) 1 - /3j(-^^xo (^f34'\f3J^''>m + (3h 



If the integral in Eq. (|132p converges for any N , by per- 
forming saddle point integration we see that the saddle 
point m'^P is solution of the equation 

m"P = mo (^(34^\(3J^^'>m'P + /3ft) , (134) 

so that, if the stability condition 

l-/3j(^)xo (/3J^^\/3j(^)m^P-K/3ft) >0, 

is satisfied, in the thermodynamic limit we arrive at the 
following expression for the free energy density // of the 
related Ising model 



^J^''>m^ + f3fo[f34'\pJ^'^m + ph) 



(135) 



Similarly, any correlation function Cj of the related Ising 
model is given in terms of the correlation function C'o of 
the unperturbed model by the following relation 



Cj = Co (/34^^ , pj^'^m + /3ft) 



(136) 
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Of course, the saddle point solution rrf^ represents 
the magnetization of the related Ising model, as can be 
checked directly by deriving Eq. (|135p with respect to 
f3h and by using Eq. PM]) . 

If the saddle point equation (|134[) has more stable solu- 
tions, the "true" free energy and the "true" observable of 
the related Ising model will be given by Eqs. (|135p and 
(I136p . respectively, calculated at the saddle point solu- 
tion which minimizes Eq. (fT35)) itself and that we will 
indicate with mj. 

Let us call /3^q'' the inverse critical temperature of the 
unperturbed model with coupling Jq^'' and zero external 
field, possibly with /J'q'' = cx3 if no phase transition exists. 
As stressed in Sec. IIIB, for the unperturbed model we 
use the expression "critical temperature" for any temper- 
ature where the magnetization mo at zero external field 
passes from to a non zero value, continuously or not. 
Note that, as a consequence, if J^-* < 0, we have for- 
mally P^Q^ — oo, independently from the fact that some 
antiferromagnetic order may be not zero. 

Let us start to make the obvious observation that a 
necessary condition for the related Ising model to have a 
phase transition at h = and for a finite temperature, 
is the existence of some paramagnetic region P/ where 
nil = 0. We see from the saddle point equation (|134p 
that, for /i = 0, a necessary condition for m/ = to 
be a solution is that he (3 < [3\^q for any (3 in P/, from 
which we get also /3c < /3^q'' . In a few lines we will see 
however that the inequality must be strict if /3^q' is finite, 
which in particular excludes the case Jq < (for which 
the inequality to be proved is trivial). 

Let us suppose for the moment that be /3c < /Jco"*- 
For (3 < /3^g'' and h = 0, the saddle point equation (|134p 
has always the trivial solution m/ = which, according 
to the stability condition, is also a stable solution if 

1-PJ^'^Xo{p4'\o) >0- (137) 

The solution mj = starts to be unstable when 

1-PJ^'^Xo{p4'\o)^0- (138) 

Eq. (|138p , together with the constrain /Jc^' < /J^q' , gives 

the critical temperature of the related Ising model /3l^'' . 
In the region of temperatures where Eq. (|137p is vio- 
lated, Eq. (|134p gives two symmetrical stable solutions 
±m/ ^ 0. From Eq. (|138p we see also that the case 
/3c^^ = /3^q' is impossible unless be J^^^ = 0, since the 



susceptibility Xo{(3Jq\o) must diverge at /3^Q^ We have 



(I) 



therefore proved that pi^'' < /3^o^ Note that for J^"' >0 
and (3 < I3^j^^ Eq. is violated only for /3 > /3^^', 

whereas for Jq'^'' < Eq. (|137p in general may be vio- 
lated also in finite regions of the (3 axis. 

The critical behavior of the related Ising model can be 
studied by expanding Eq. (|134p for small fields. How- 
ever, we find it more convenient to expand L(m) in series 



around to = since in this way everything can be cast 
in the standard formalism of the Landau theory of phase 
transitions. From Eq. (|133p . taking into account that 
the function xo iPJojPh) is an even function of (3h, we 
have the following general expression valid for any to, /3 
and small h 

L{m) = 13 h (fi4'\o) - mo (/3J^'\ O) /3/i + V (mX;39) 

where we have introduced the Landau free energy density 
'ip{m) given by 



Tp (to) 



1 



1. 



-am 



-bm 



1 



-cm 



where 



1 - 



nPh + A (/3/o) (/3jy\ /3J(^)to) , (140) 
f3J^'^Xo{p4'\o)]f3J^'\ (141) 



92 



9(/3/i)2 
' d{l3h)^ 



Xo 



3! 



5! 



(142) 



(143) 



h^mo [f34'\o) J^'^+Xo {p4'\o) /3(^)j(^)/3/i(144) 

finally, the last term A(/3/o) (^(3,4'\ (3J^^^m^ is defined 

implicitly to render Eqs. p39p and p40p exact, but 
terms 0{h^) and 0{m^h); explicitly 



A 



{I3h)(f34'\pj^'^ 



E 

fc=4 



d{(3h) 



2k~2 



Xo 



2k 



■fl45) 



Finally, to come back to the original random model, 
we have just to perform the reversed mapping transfor- 
mations (fT29)) and p30)) in Eqs. p33)) - (pl5l) . As a result 
we get immediately Eqs. (Il5])-(|13|), but Eq. (HH). 



IX. DERIVATION OF EQ. ([21]) AND EQS. 

HMD-dZSI) 

Concerning Eq. (PT|) for the full expression of the free 
energy density, it can be obtained by using Eqs. (jlOSp . 
(frog)) , and (fTTT)) . Here Lpi is the high temperature 

part of the free energy density of the related Ising model 
we have just solved: 



1 



I3fi = lim — V log cosh(/3J($ 

N^OQ I\ ^ — ^ 



N -1 



(ij)6ro 



X log 



cosh 



(^/3j(-fV^)] +log[2cosh(/3/i)] +<^7(146) 
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where we have taken into account the fact that our re- 
lated Ising model has IFqI connections with coupling Jg^^ 
and N(N — l)/2 connections with the coupling J'^^^/A^. 
By using Eq. (|135p calculated in rrij and Eq. (|146p . for 
large N we get 



lim — 

N->oo N ^ 

(i,i)ero 



log 



cosh(/3jy^) 



log [2cosh(/3/i)] +0 



(147) 



limit we have (/Sq (/3 J'o , 0) = 0; if instead 2 < do < cxd, 
93o(/3'^o, 0) ^ 0. We see here therefore what anticipated 
in Sec. VIIA: when Jq 7^ 0, the symmetry among the 
random couplings is broken and for do sufficiently high 
this reflects in a non zero (p*-^' also in the P region. 

Next we prove Eqs. ([Ml) - ([7n)) . To this aim we have 
to calculate Eq. (|150p at the leading solution mS^^ and 
to compare 93^^) and (/j^^*^'. Note that the term in the 
square parenthesis of Eq. (|150|) is non negative since 
rrS^^ is the absolute minimum of L(^). We recaU that 
for critical temperature we mean here any temperature 
lying on the boundary P-F or P-SG, so that 
for any /3 in the P region. 



Therefore, on using Eq. (jll5p . for the non trivial part 
(/?(^) of the random system, up to corrections 0(l/iV), 
we arrive at 

^(S) ^ (m(^))' -log[2cosh(/3;i)] 



lim — 

(i,i)Gro 



log 



cosh(/34^^ 



-/3/o(/?jf\/3J(''^m(^)+/3/i) (148) 
In terms of the function L(^^(to) Eq. p48p reads as 



1 



lim 

N^oa N 



J2 log [cosh(/3jf ^ 
(ij)ero 
- log [2 cosh(/3/?,)] . 



(149) 



By using Eqs. (fTOS]) . (fMO]) . (fTT5ll and ITTTI) . with I = 1 
or 2 for E=F or E=SG, respectively, we get Eq. (PT|) . 
For /i = Eq. (|149p can conveniently be rewritten also 

as 



(/3jf\0) 
L(^) (0) - L(^) (m(^; 



(150) 



where 



^o(/?Jo,/?M = -/9/o(/3jf\/3/i 

log \cosh{(3jP) 



1 



— lim — > 
-log [2cosh(/3/i)] , 



(151) 



is the high temperature part of the free energy density 
of the unperturbed model with coupling Jq and exter- 
nal field h. There are some important properties for the 
function (Pq^PJq, 0): it is a monotonic increasing function 
of pJo; if the lattice £0 has only loops of even length, 
ipo{f3Jo,0) is an even function of /3Jo; furthermore, if 
do < 2 and the coupling-raiige is finite, or if do = 00 
at least in a wide sense jlS], in the thermodynamic 



A. Jo > 

If Jo > 0, for both the solution with label F and SG, 
we have only one second-order phase transition so that 
771(F) = and jfi'^^) = 0, respectively, are the stable 
and leading solutions even on the boundary with the P 
region. 

Let be (po{-,0) ^ 0. From 



T ^- /3(F) ^ a(SG) 

Let us suppose pc < Pc ■ 



Eq. pSOp and by using Jq > Jq , we see that 



(F)| 



Vo 



(F) t(SG) 



,0 



(152) 



Finally, by using this result and the general rule given 
by Eqs. (|115p and (|117p . we see (and with a stronger 
reason, due to the factor 1/2 appearing in these equations 
for the SG solution) that the stable phase transition is 
the P-F one: /3c = P', 
for < f3 < (3, 



Similarly, by using Eq. pSO]) 
we see that even for any (3 in 



the interval (/3c , /3c ) the stable solution is that with 
label F. This last observation makes also clear that if 
</3o(', 0) = we reach the same conclusion: F is the stable 

/pN (SCI 

phase in all the region /3c < (3 < (3c and in particular 
this implies also that the stable phase transition is the 
P-F one: (3c = (3^ . 

Let us suppose (3c > (3c . If v3o('j 0) ^ 0, we arrive 

at 



/3(SG)j(SG) 







(153) 



Finally, by using this result and the general rule given 
by Eqs. (jllSp and (|117p . we see (and with a stronger 
reason) that the stable phase on the boundary is that 
predicted by the F solution which has zero magnetization 



at P, 



(SG) 



This does not imply that (3c 



Pc^\ but only 



that /3f < (3c< /3r^ If instead V9o(-,0) = 0, by using 
Eq. dHOD for < /3 < (3P, we see that for any 

f3 in the interval {(3c ,/5c ) the stable solution is SG 
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and then, in particular, the stable boundary is P-SG: 

/3c = /3f°^ 



B. Jo < 

If Jo < 0, for the solution with label F, we may have 
both first and second-order phase transitions. In the first 
case we cannot in general assume that is the stable 
and leading solution on the boundary with the P region: 
i^^^^l/^m ^ in general. As a consequence, for a first- 
order transition the term in square parenthesis of Eq. 
()150p may be non zero even on the critical surface. Fur- 
thermore, as Jo < 0, for the solution F we have at least 
two critical temperatures that we order as > p'^^ ■ 
However, despite of these complications, if we assume 
that Co has only loops of even length, (po{-,0) turns 
out to be an even function and, due to the inequality 
I Jg^^ I > Jq^'^'' , almost nothing changes in the arguments 
we have used in the previous case Jo > 0. 

Let us consider first the surfaces and /Jc^*^^ . Inde- 
pendently of the kind of phase transition, first or second- 
order, we arrive again at Eqs. (|152p and (|153p . for /3^2'' < 
/3c^^'' and (3'^^ > pi^'^^ , respectively, with the same pre- 
scription for the cases (po{-,0) ^ 0, or ipo{-,0) = 0. 

Let us now consider the surfaces /3^P and f3i^'^\ If 
P^-^^' < Pc and </Jo(-,0) ^ 0, for any /3 in the interval 
[/?cf\/3c^°^] we have, 

^(^)|, = ^o (/3jr,o) 

>^(^^^\,^^o{pjr\o), (154) 

SO that the interval [P^^ ,Pc ] is a stable P region corre- 
sponding to the solution with label F. Similarly, we arrive 
at the same conclusion if (/3o(',0) = 0. However, the in- 
terval of temperatures where the P region [P^^\pf'^^] 
is stable can be larger when > 2. In fact (exactly 
as we have seen for Jo > 0) in this case the P-SG stable 
boundary may stay at lower temperatures. Finally, let us 
consider the case /J^^ > Pc ■ If i^sol', 0) = 0, by using 
Eq. (|150p we see that for any P > /SF*^' we have that 
the stable solution corresponds to the SG one, so that 
there is no stable boundary with the P region. If instead 
ipo{-,0) Oj due to the fact that (p^^^^lfs and ^^^^1/3 grow 
in a different way with /3, we are not able to make an ex- 
act comparison, and it is possible that the P-F boundary 

(F) 

becomes stable starting from some Pd with Pd > P^i ■ 
In general, as in the case Jq < 0, we could have one (or 
even more) sectors where the P region corresponding to 
the solution with label F is stable. 



X. CONCLUSIONS 

In this paper we have presented a novel and general 
method to analytically face random Ising models de- 
fined over small- world networks. The key point of our 
method is the fact that, at least in the P region, any 
such a model can be exactly mapped to a suitable fully 
connected model, whose resolvability is in general non 
trivial for do > 1, but still as feasible as a non random 
model. As a main result we then derive a general self- 
consistent equation, Eq. ([T3]), which allows to describe 
effectively the model once the magnetization of the un- 
perturbed model in the presence of a uniform external 
field, 'mo{pJo, ph), is known. 

The physical interpretation of this general result is 
straightforward. From Eqs. (jlSp we see that, concerning 
the magnetization m^^\ the effect of adding long range 
Poisson distributed bonds implies that the system - now 
perturbed - feels, besides the coupling Jo, also an effec- 
tive external field J^^^ shrunk by m^^^ itself. Concerning 
m^^'^\ the effect is that the system now feels a modified 
effective coupling Jq^'^'' and an effective external field 
j(SG) shrunk by m'^^'^'> itself. 

We are therefore in the presence of an effective field the- 
ory which, as opposed to a simpler mean-field theory, de- 
scribes mS^^ and m*^^*^^ in terms of, not only an effective 
external field, but also through the non trivial function 
nT-oiPJoTPh) which, in turn, takes into account the cor- 
relations due to the non zero short-range coupling Jo or 
Jq felt by the unperturbed system. The combination 
of these two effects gives rise to the typical behavior of 
models defined over small- world networks: the presence 
of a non zero effective external field causes the existence 
of a phase transition also at low do dimension. However, 
the precise determination of both the critical surface and 
the correlation functions is obtained in a non trivial way 
via the unperturbed magnetization mo{pJoT Ph). 

We have used the method introduced to analyze the 
critical behavior of generic models with Jo > and 
Jo < and showed that they give rise to two strictly 
different phase transition scenarios. In the first case, we 
have a mean-field second-order phase transition with a 
finite correlation length. Whereas, in the second case, 
we obtain multiple first- and second-order phase tran- 
sitions. Furthermore, we have shown that the com- 
bination of the F and SG solutions results in a total 
of four possible kinds of phase diagrams according to 
the cases (1) (Jo > 0; do < 2, or do = oo), (2) 
(Jo > 0; 2 < do < oo), (3) (Jo < 0; do < 2, or do = oo), 
and (4) (Jo < 0; 2 < do < oo). One remarkable dif- 
ference between systems with do < 2, or do = oo and 
those with 2 < do < oo, is that in the latter case we 
have, in principle, also first-order P-SG phase transitions 
and, moreover, re-entrance phenomena are in principle 
possible even for Jo > 0. 

In Sees. IV, V, and VI we have applied the method 
to solve analytically those models for which the unper- 
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turbed magnetization mo(/3Jo, f3h) is known analytically, 
i.e., the small-world models in dimension dp — 0,1, oo, 
corresponding to an ensemble of non interacting units 
(spins, dimers, etc.), the one dimensional chain, and 
the spherical model, respectively. In particular, we have 
studied in detail the small-world model defined over the 
one dimensional chain with positive and negative short- 
range couplings showing explicitly how, in the second 
case, multicritical points with first- and second-order 
phase transition arise. Finally, the small-world spheri- 
cal model - an exact solvable model (in our approach) 
with continuous spin variables - has provided us with 
an interesting case study to explore what happens as do 
changes continuously fi:om to oo. As expected from 
general grounds, unlike the non random version of the 
model, the small-world model presents always a finite 
temperature phase transition, even in the limit do — > 0^. 
This latter result - besides to be consistent with what 
we have found in the do = 0-dimensional discrete mod- 
els - has a simple physical explanation in our approach. 
In fact, it consists on mapping the small- world model (a 
random model) to a corresponding non random model 
(no long-range bonds) but immersed in an effective uni- 
form external field which is active as soon as the added 
random connectivity c is not zero (see Eqs. (fT5|) -(fT7 | ). 

Many interesting variants of the above models can be 
considered and are still analytically solvable by our ap- 
proach (see also the generalizations considered in Sec. 
HIE). However, our approach can be also applied nu- 
merically to study more complex small-world models for 
which the corresponding unperturbed model is not ana- 
lytically available [29|. In fact, the numerical complexity 
in solving such small-world models is comparable to that 
required in solving a non random model immersed in a 
uniform external field. 

Models defined on complex small-world networks [s^l 
are an interesting subject of future work. 
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couplings Jo but now immersed in an arbitrary (non ho- 
mogeneous) external field {/in}, where n — 1, . . . , N . Af- 
ter using the Gaussian transformation we have the fol- 
lowing partition function: 



(Al) 



where we have introduced the function 

L{m) = ^ Jm2 + 13 fo (/3Jo, {PJm + PK}) , (A2) 

/o iPJo, {Phn}) being the free energy density of the un- 
perturbed model in the presence of an arbitrary external 
field {(3hn}. By using 



dph, /3/o(/3Jo, {Ph^}) - -mo^ipJo, {/3/i„}), (A3) 



and 



def 

XO;i.]{PJo,{f^h,i}) = {(Ji<Jj)o - {<Ji)o{crj)o 



di3h, mo.(/3Jo,{/3/i„})(A4) 



we get 
L'{m) = PJ 



CA5) 



L"{m) = p.J 



1-PJ 



N 



X ^Xo^j {PJo,{PJ'm' + Pl^n}) 



'1,3 



(A6) 



By performing the saddle point integration we see that 
the saddle point rrf^ is solution of the equation 

m^P = l^mo,(/3Jo,{/3Jm^P + /3/i„}) (A7) 
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hence, by using 

Xo {PJoAPJm + pK}) = 



(A8) 



APPENDIX A: GENERALIZATION TO NON 
HOMOGENEOUS EXTERNAL FIELD 



In this appendix we prove Eq. ([55]) calculating the 
0{^/N) correction responsible for the divergence of the 
susceptibility of the random system at Tc. To this aim we 
firstly need to generalize our method to an arbitrary ex- 
ternal field. Let us consider again a fully connected model 
having - as done in Sec. VIII - long-range couplings J 
(for brevity we will here omit the label /) and short-range 



we see that if the stability condition 

1 - pjxo iPJoAPJ^"^ + Phn}) > 0, 



(A9) 



is satisfied, in the thermodynamic limit we arrive at the 
following expression for the free energy density / of the 
related Ising model immersed in an arbitrary external 
field 



^ Jm^ + Pfo iPJo, {pJm + PK}) 



(AlO) 
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On the other hand, by derivation with respect to (3hi and 
by using Eq. (jA7p . it is immediate to verify that 

m, ''^ (a,) = mo,(/3Jo, {/3Jm^P + PK}), (All) 

and then also (from now on for brevity on we omit the 
symbol ^'^) 



(A12) 



We want now to calculate the correlation functions. 
From Eq. (jAll[) . by deriving with respect to phj we 
have 

i^mu:\ = ^Xo;i,i{PJo,{PJ'm + PK}) 



d{l3hj 
f ^ dm ^ 



(A13) 



which, by summing over the index i and using (jA12p . 
gives 



9m j^Y.iXo;i:]{PJo,{PJ'm + /3hn}) 



(A14) 



We can now insert Eq. (|A14p in the rhs of Eq. (IA13|) to 

get 



Xij = XO;i,jiPJo,{PJ'm + Phn}) 

^ N 1- pJxoiPJo. {PJm + PK}) ' 



(A15) 



where for brevity we have omitted the argument in Xo;i,j 
and xo;i,fcj which is the same of xo appearing in the de- 
nominator. If we now come back to choice a uniform 
external field hn = h, n = 1, . . . , we can use trans- 
lational invariance and for the related Ising model (fully 
connected) we obtain the following correlation function 



PJ [xo{pJo,pjm + ph)f 
N 1- pJxo{pJa,PJm + Ph) 
XO;ij{PJo,PJ'm- + Ph). 



(A16) 



Finally, by performing the mapping substitutions (|129p 
and (I130|) we arrive at Eq. ([33|) . 

Similarly, any correlation function C of the related 
Ising model will be given by a similar formula with the 
leading term Cq plus a correction 0(1/A^) becoming im- 
portant only near Te- 
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